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METHODS FOR THE COMPENSATION OF IMAGING TECHNIQUE IN THE 
PROCESSING OF RADIOGRAPHIC IMAGES 

Field of the Invention 

[0001] The present invention is in the field of imaging and analysis thereof. In 
5 particular, methods and compositions for accurately analyzing images to determine 
bone mineral density and/or bone structure are described. 

Background of the Invention 

[0002] Osteoporosis is a condition that affects millions of Americans. 
Osteoporosis refers to a condition characterized by low bone mass and 
10 microarchitectural deterioration of bone tissue, with a consequent increase of bone 
fragility and susceptibility to fracture. Osteoporosis presents commonly with vertebral 
fractures or hip fractures due to the decrease in bone mineral density and deterioration 
of structural properties and microarchitecture of bone. 

[0003] Imaging techniques are important diagnostic tools, particularly for bone 
15 related conditions. Currently available techniques for the noninvasive assessment of 
the skeleton for the diagnosis of osteoporosis or the evaluation of an increased risk of 
fracture include dual x-ray absorptiometry (DXA) (Eastell et al. (1998) New Engl J. Med 
338:736-746); quantitative computed tomography (QCT) (Cann (1988) Radiology 
166:509-522); peripheral DXA (pDXA) (Patel et al. (1999) J Clin Densitom 2:397-401); 
20 peripheral QCT (pQCT) (Gluer et. al. (1997) Semin Nucl Med 27:229-247); x-ray image 
absorptiometry (RA) (Gluer et. al. (1997) Semin Nucl Med 27:229-247); and 
quantitative ultrasound (QUS) (Njeh et al. "Quantitative Ultrasound: Assessment of 
Osteoporosis and Bone Status" 1999, Martin-Dunitz, London England; U.S. Patent No. 
6,077,224, incorporated herein by reference in its entirety). (See, also, WO 9945845; 
25 WO 99/08597; and U.S. Patent No. 6,246,745). 

[0004] DXA of the spine and hip has established itself as the most widely used 
method of measuring BMD. Tothill, P. and D.W. Pye, (1992) Br J Radiol 65:807-81 3. 
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The fundamental principle behind DXA is the measurement of the transmission through 
the body of x-rays of 2 different photon energy levels. Because of the dependence of 
the attenuation coefficient on the atomic number and photon energy, measurement of 
the transmission factors at 2 energy levels enables the area densities (i.e., the mass per 

5 unit projected area) of 2 different types of tissue to be inferred. In DXA scans, these are 
taken to be bone mineral (hydroxyapatite) and soft tissue, respectively. However, it is 
widely recognized that the accuracy of DXA scans is limited by the variable composition 
of soft tissue. Because of its higher hydrogen content, the attenuation coefficient of fat is 
different from that of lean tissue. Differences in the soft tissue composition in the path 

10 of the x-ray beam through bone compared with the adjacent soft tissue reference area 
cause errors in the BMD measurements, according to the results of several studies. 
Tothill, P. and D.W. Pye, (1992) Br J Radiol, 65:807-813; Svendsen, O.L., et al., (1995) 
J Bone Min Res 10:868-873. Moreover, DXA systems are large and expensive, ranging 
in price between $75,000 and $150,000. 

15 [0005] Quantitative computed tomography (QCT) is usually applied to measure 
the trabecular bone in the vertebral bodies. Cann (1988) Radiology 166:509-522. QCT 
studies are generally performed using a single kV setting (single-energy QCT), when 
the principal source of error is the variable composition of the bone marrow. However, 
a dual-kV scan (dual-energy QCT) is also possible. This reduces the accuracy errors 

20 but at the price of poorer precision and higher radiation dose. Like DXA, however, 
QCT are very expensive and the use of such equipment is currently limited to few 
research centers. 

[0006] Quantitative ultrasound (QUS) is a technique for measuring the peripheral 
skeleton. Njeh et al. (1997) Osteoporosis Int 7:7-22; Njeh et al. Quantitative Ultrasound: 
25 Assessment of Osteoporosis and Bone Status. 1999, London, England: Martin Dunitz. 
There is a wide variety of equipment available, with most devices using the heel as the 
measurement site. A sonographic pulse passing through bone is strongly attenuated as 
the signal is scattered and absorbed by trabeculae. Attenuation increases linearly with 
frequency, and the slope of the relationship is referred to as broadband ultrasonic 
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attenuation (BUA; units: dB/MHz). BUA is reduced in patients with osteoporosis 
because there are fewer trabeculae in the calcaneus to attenuate the signal. In addition 
to BUA, most QUS systems also measure the speed of sound (SOS) in the heel by 
dividing the distance between the sonographic transducers by the propagation time 

5 (units: m/s). SOS values are reduced in patients with osteoporosis because with the 
loss of mineralized bone, the elastic modulus of the bone is decreased. There remain, 
however, several limitations to QUS measurements. The success of QUS in predicting 
fracture risk in younger patients remains uncertain. Another difficulty with QUS 
measurements is that they are not readily encompassed within the WHO definitions of 

10 osteoporosis and osteopenia. Moreover, no intervention thresholds have been 
developed. Thus, measurements cannot be used for therapeutic decision-making. 
[0007] There are also several technical limitations to QUS. Many devices use a 
foot support that positions the patient's heel between fixed transducers. Thus, the 
measurement site is not readily adapted to different sizes and shapes of the calcaneus, 

15 and the exact anatomic site of the measurement varies from patient to patient. It is 
generally agreed that the relatively poor precision of QUS measurements makes most 
devices unsuitable for monitoring patients' response to treatment. Gluer (1997) J Bone 
Min Res 12:1280-1288. 

[0008] Radiographic absorptiometry (RA) is a technique that was developed 
20 many years ago for assessing bone density in the hand, but the technique has recently 
attracted renewed interest. Gluer et al. (1997) Semin Nucl Med 27:229-247. With this 
technique, BMD is measured in the phalanges. The principal disadvantage of RA of the 
hand is the relative lack of high turnover trabecular bone. For this reason, RA of the 
hand has limited sensitivity in detecting osteoporosis and is not very useful for 
25 monitoring therapy-induced changes. 

[0009] Peripheral x-ray absorptiometry methods such as those described above are 
substantially cheaper than DXA and QCTwith system prices ranging between $15,000 
and $35,000. However, epidemiologic studies have shown that the discriminatory ability 
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of peripheral BMD measurements to predict spine and hip fractures is lower than when 
spine and hip BMD measurements are used. Cummings et al. (1993) Lancet 341 :72- 
75; Marshall et al. (1996) Br Med J 312:1254-1259. The main reason for this is the lack 
of trabecular bone at the measurement sites used with these techniques. In addition, 

5 changes in forearm or hand BMD in response to hormone replacement therapy, 
bisphosphonates, and selective estrogen receptor modulators are relatively small, 
making such measurements less suitable than measurements of principally trabecular 
bone for monitoring response to treatment. Faulkner (1998) J Clin Densitom 1:279- 
285; Hoskings et al. (1998) N Engl J Med 338:485-492. Although attempts to obtain 

10 information on bone mineral density from dental x-rays have been attempted (See, e.g., 
Shrout et al. (2000) J. Periodonol. 71:335-340; Verhoeven et al. (1998) Clin Oral 
Implants Res 9(5):333-342), these have not provided accurate and reliable results. 
[0010] Furthermore, current methods and devices do not generally take into account 
bone structure analyses. See, e.g., Ruttimann et al. (1992) Oral Surg Oral Med Oral 

15 Pathol 74:98-1 10; Southard & Southard (1992) Oral Surg Oral Med Oral Pathol 73:751- 
9; White & Rudolph, (1999) Oral Surg Oral Med Oral Pathol Oral Radiol Endod 88:628- 
35. 

[001 1] Thus, although a number of devices and methods exist for evaluating 
bone, there are a number of limitations on such devices and methods. Consequently, 

20 the inventors have recognized the need, among other things, to provide methods and 
compositions that result in the ability to obtain accurate bone mineral density and bone 
structure information from images (e.g., radiographic images) and data. 

Summary of the Invention 
[0012] In one aspect, the disclosure provides a method to derive information 

25 regarding one or more bone parameters from an image, the method comprising the 
steps of: (a) obtaining an image comprising bone from a subject; (b) defining two or 
more regions of interest (ROIs) in the image; and (c) analyzing a plurality of positions in 
the ROIs to determine one or more parameters selected from the group consisting of 
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bone microarchitecture, bone macro-anatomy, biomechanical parameters and 
combinations thereof of the ROIs. In certain embodiments, the ROIs are overlapping. 
The positions analyzed in the ROIs may be at regular intervals relative to one another 
or, alternatively, may be irregularly spaced relative to each other. Thus, in certain 
5 embodiments, the methods involve determining bone micro-architecture, for example by 
analyzing positions at regular intervals. In other embodiments, the methods involve 
determining bone macro-anatomy, for example by analyzing positions at irregular 
intervals in the image. 

[0013] In any of the methods described herein, the image can be two- 
10 dimensional (2D) or three-dimensional (3D). The images may be x-rays, MRI images, 
CAT scan images, or any other image including bone. In any of the methods, the 
image may be an electronic image. 

[0014] In any of the methods described herein, the subject can be, for example, 
an osteoporosis subject. 

15 [0015] In another aspect, this disclosure relates to a method of generating a map 
of one or more bone parameters, the method comprising the steps of (a) obtaining 
information on bone parameters according to the method of any of methods described 
herein; and (b) identifying regions of the image that exhibit similar parameter 
characteristics, thereby creating a parameter map of the image. 

20 [0016] In yet another aspect, a method of predicting a fracture path in a subject is 
provided, the method comprising the steps of: (a) generating multiple parameter maps 
according to any of the methods of generating parameters maps described herein; (b) 
generating a composite parameter map from the multiple parameters maps of step (a); 
and (c) analyzing the composite parameter map to identify possible fracture paths. 

25 [0017] In yet another aspect, the invention includes a method of predicting a 
fracture path in a subject, the method comprising the steps of: (a) analyzing of one or 
more parameter maps preparing according to any of the methods described herein, 
wherein the analysis is watershed segmentation analysis or Markov random field 
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analysis; and (c) identifying possible fracture paths based on the analysis of step (a), 
thereby predicting a fracture path in the subject. 

[0018] In another aspect, the invention includes a method of predicting the risk of 
fracture in a subject, the method comprising the steps of: (a) generating a finite element 
5 model from one or more parameter maps obtained according any of the methods 

described herein; (b) applying simulated force vectors that would occur during a fracture 
incident to the finite element model generated in step(s); and (c) determining the 
minimum forces required for fracture to occur, thereby estimating the risk of fracture. 

[0019] In a still further aspect, the invention includes a method of determining the 
10 risk of fracture in a subject comprising: (a) predicting a fracture path according to any of 
the methods of predicting fracture path as described herein; (b) evaluating one or more 
selected bone parameters along the predicted fracture path, thereby estimating the risk 
of fracture. 

[0020] In another aspect, the invention includes a method of treating a subject 
15 with bone disease comprising (a) obtaining an image from a subject; (b) analyzing the 
image obtained in step (a) using any of the methods described herein; (c) diagnosing a 
bone disease based on the analysis of step (b); and (d) selecting and administering a 
suitable treatment to said subject based on said diagnosis. 

Brief Description of the Drawings 

20 [0021] Fig. 1 shows an example of a dental x-ray. A calibration phantom 1 10 is 
seen. Regions of interest 120 have been placed for measurement of bone mineral 
density or structure. 

[0022] Fig. 2 shows another example of a dental x-ray. A calibration phantom 
1 10 is seen. Regions of interest 120 have been placed for measurement of bone 
25 mineral density or structure. 

[0023] Fig. 3 shows an example of an analysis report resulting from a 
measurement of mandibular or maxillary bone mineral density. A subject (X) is more 
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than one standard deviation below the mean of age-matched controls (x-axis age, y- 
axis arbitrary units BMD). 

[0024] Fig. 4 shows an example of a V-shaped calibration phantom 1 1 0 mounted 

on a tooth 120. Gums are also shown 130. 
5 [0025] Fig. 5 shows an example of a holder 115 for a calibration phantom 1 1 0. 

The holder 115 is mounted on a tooth 120. Gums are also shown 130. 

[0028] Fig. 6, panels B through E shows gray value profiles along different rows 

of pixels used for locating dental apices. From top to bottom, the characteristic peaks 

for the dental roots (shown in dental x-ray panel A) gradually disappear. 
10 [0027] Fig. 7 shows a Hough transform (panel A) of a test image (panel B). All 

collinear points from the same line are transformed into sinusoidal curves that intersect 

in a single point (circles). 

[0028] Fig. 8 shows a Hough transform (panel A) of a skeletonized trabecular 

bone x-ray image (panel B). The white regions in panel A indicate longer segments and 
15 predominant angles. 

[0029] Fig. 9 shows the effect of varying size of structuring element E 2 ; 

calibration phantom image with lines of varying width (1, 3, 5, 7, 9, 11, 13 pix) (top left); 

skeleton operation performed using E 2 with a diameter of 3 pix (top right), 7 pix (bottom 

left), and 11 pix (bottom right), respectively. 
20 [0030] Fig. 10 shows the effect of varying size of structuring element E 2 ; gray 

scale image of trabecular bone (top left, panel A); skeleton operation performed using 

E 2 with a diameter of 3 pix (top right, panel B); 7 pix (bottom left, panel C) and 1 1 pix 

(bottom right, panel D), respectively. 

[0031] Fig. 1 1 shows gray value surface plot of an anatomical region of interest 
25 from a dental x-ray (inset) used for fractal analysis. 

[0032] Fig. 1 2 shows an example of a hygienic cover holder that includes 
compartments for a calibration phantom and a fluid-filled bolus back. 
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[0033] Fig. 1 3 shows an example of an anatomical region of interest (black dot), 
determined relative to the teeth or to the convexity/concavity of the mandible. 

[0034] Fig. 14 shows an example of three snatomical region of interests (black 
dots), determined relative to the teeth or to the convexity/concavity of the mandible. 

5 [0035) Fig. 15 is a side view of an exemplary system for minimizing tube 

angulation as described herein. In the Figure, the system is shown as a dental x-ray 
system. An extension tubing (200) is attached to a ring-shaped Rinn holder (102). The 
outer diameter of the extension tubing is slightly smaller than the inner diameter of the 
tube located in front of the dental x-ray system/dental x-ray tube. The extension tubing 
10 can then be inserted into the metal tube thereby reducing tube angulation and resultant 
errors in bone apparent density and bone structural measurements. 

[0036] Fig. 16 depicts an example of a regular interval sampling field for 
microarchitecture (+) and a higher density sampling field for macro-anatomical 
features^) on a femur radiograph. White rectangles are examples of overlapping 
15 window positioning. 

[0037] Fig. 17 depicts watershed segmentation boundaries superimposed on a 
parameter map. The two white lines are the actual fracture paths resulted from an in- 
vitro mechanical loading test. 

[0038] Fig. 18 is a flowchart depicting an exemplary process to determine 
20 fracture risk using overlapping window processing and fracture paths prediction. 

[0039] Fig. 19 depicts a Markov random field analysis by modeling particular joint 
feature distributions as they are estimated at each image element or image 
neighborhood. 

[0040] Fig. 20 depicts an exemplary model definition for trabecular pattern 
25 density characterization in a region of interest (ROI) with a noise model P(N) and 
characteristic structure pattern given a density level P(l | Ti). 

[0041] Fig. 21 depicts exemplary Bayes' Rule analysis. 
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[0042] Fig. 22 depicts an example of a regular interval sampling field for 
microarchitecture (+) and a higher density sampling field for macro-anatomical features 
(*) on a spine radiograph. White rectangles are examples of overlapping window 
positioning. 

5 [0043] Fig. 23 depicts an example of a sampling field of varying density for 
microarchitecture (+, x, diamond) and a regular sampling field for macro-anatomical 
featuresO on a knee radiograph. White rectangles are examples of overlapping window 
positioning. 

[0044] Fig. 24 depicts an example of an application of structure extraction and 
10 measurement for therapeutic monitoring using spine x-ray. White outline of extracted 
structure are show in (a) before treatment, and (b) after treatment. 

Detailed Description of the Invention 

[0045] The following description is presented to enable any person skilled in the 
art to make and use the invention. Various modifications to the embodiments described 

15 will be readily apparent to those skilled in the art, and the generic principles defined 

herein can be applied to other embodiments and applications without departing from the 
spirit and scope of the present invention . Thus, the present invention is not intended to 
be limited to the embodiments shown, but is to be accorded the widest scope consistent 
with the principles and features disclosed and shown herein. To the extent necessary 

20 to achieve a complete understanding of the invention disclosed, the specification and 
drawings of all issued patents, patent publications, and patent applications cited in this 
application are incorporated herein by reference. 

[0046] The practice of the present invention employs, unless otherwise indicated, 
conventional methods of imaging and image processing within the skill of the art. Such 
25 techniques are explained fully in the literature. See, e.g., WO 02/22014, incorporated 
herein in its entirety by reference; X-Ray Structure Determination: A Practical Guide, 
2nd Edition, editors Stout and Jensen, 1989, John Wiley & Sons, publisher; Body CT: A 
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Practical Approach, editor Slone, 1999, McGraw-Hill publisher; The Essential Physics of 
Medical Imaging, editors Bushberg, Seibert, Leidholdt Jr & Boone, 2002, Lippincott, 
Williams & Wilkins; X-ray Diagnosis: A Physician's Approach, editor Lam, 1998 
Springer-Verlag, publisher; and Dental Radiology: Understanding the X-Ray Image, 

5 editor Laetitia Brocklebank 1 997, Oxford University Press publisher. 

[0047] Methods and compositions useful in analyzing images are described. In 

particular, the invention includes methods of obtaining and/or deriving information about 
bone mineral density and/or bone structure from an image. Additionally, the present 
invention relates to the provision of accurate calibration phantoms for use in determining 

10 bone structure and methods of using these calibration phantoms. In particular, the 
present invention recognizes for the first time that errors arising from misplacement of 
interrogation sites in dental or hip x-rays of bone density and/or bone structure can be 
corrected by positioning the x-ray tube, the detector and/or the calibration reference with 
respect to an anatomical landmark (or anatomical region of interest). 

15 [0048] Advantages of the present invention include, but are not limited to, (i) 
providing accessible and reliable means for analyzing x-rays; (ii) providing non-invasive 
measurements of bone structure and architecture and macro-anatomy; (iii) providing 
methods of diagnosing bone conditions {e.g., osteoporosis, fracture risk); (iv) providing 
methods of treating bone conditions; and (iv) providing these methods in cost-effective 

20 manner. 

1.0. Obtaining Data from Images 
[0049] An image can be acquired using well-known techniques from any local 
site. Non-limiting examples of imaging techniques suitable for acquiring images from 
which data can be obtained include, ultrasound, CAT scan, MRI and the like. See, also, 
25 "Primer of Diagnostic Imaging," 3rd edition, eds. Weissleder et al. (2002), Mosby Press; 
and International Publication WO 02/22014. 

[0050] In certain aspects, 2D planar x-ray imaging techniques are used. 2D 
planar x-ray imaging is a method that generates an image by transmitting an x-ray 
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beam through a body or structure or material and by measuring the x-ray attenuation on 
the other side of said body or said structure or said material. 2D planar x-ray imaging is 
distinguishable from cross-sectional imaging techniques such as computed tomography 
or magnetic resonance imaging. If the x-ray image was captured using conventional x- 

5 ray film, the x-ray can be digitized using any suitable scanning device. Digitized x-ray 
images can be transmitted over a networked system, e.g. the Internet, into a remote 
computer or server. It will be readily apparent that x-ray images can also be acquired 
using digital acquisition techniques, e.g. using photostimulable phosphor detector 
systems or selenium or silicon detector systems, the x-ray image information is already 

10 available in digital format which can be easily transmitted over a network. In other 
embodiments, 3D images are acquired, for example, using 3D imaging techniques 
and/or by creating 3D images from 2D images. 

[0051] Any images can be used including, but not limited to, digital x-rays and 
conventional x-ray film (which can be digitized using commercially available flatbed 

15 scanners). In certain embodiments, the x-ray is of the hip region, for example 

performed using standard digital x-ray equipment (Kodak DirectView DR 9000, Kodak, 
Rochester, NY). Patients are typically positioned on an x-ray table in supine position, 
parallel to the long axis of the table, with their arms alongside their body. The subject's 
feet may be placed in neutral position with the toes pointing up or in internal rotation or 

20 may be placed in a foot holder such that the foot in a neutral position (0° rotation) or in 
any desired angle of rotation (e.g., internal or external) relative to neutral (see, also 
Example 8 below). Foot holders suitable for such purposes may include, for example, a 
base plate extending from the foot, for example, from the mid to distal thigh to the heel. 
The base plate preferably sits on the x-ray table. The patients' foot is positioned so that 

25 the posterior aspect of the heel is located on top of the base plate. The medial aspect of 
the foot is placed against a medial guide connected rigidly to the base plate at a 90° 
angle by any suitable means (e.g., straps, velcro, plastic, tape, etc.). A second, lateral 
guide attached to the base plate at a 90° angle with a sliding mechanism can then be 
moved toward the lateral aspect of the foot and be locked in position, for example when 



11 



WO 2004/086972 



PCT7US2004/00<M<>5 



it touches the lateral aspect of the foot. The use of a foot holder can help improve the 
reproducibility of measurements of bone structure parameters or macro-anatomical 
and/or biomechanical parameters. 

[0052] As will be appreciated by those of skill in the art, the patient or subject can 

5 be any warm-blooded animal. Typically patients, or subjects, are chosen from the class 
Mammalia. Thus, for example, patients, or subjects, would include humans and 
nonhuman primates such as chimpanzees and other apes and monkey species; farm 
animals such as cattle, sheep, pigs, goats and horses; domestic mammals such as 
dogs and cats; laboratory animals including rodents such as mice, rats and guinea pigs, 

10 and the like. To the extent desirable, other non-mammals can be subjected to the 
protocols described herein without departing from the scope of the invention, 
[0053] Persons of skill in the art will appreciate that macro-anatomical 
parameters generally describe the shape, size or thickness of bone and/or surrounding 
structure. Oftentimes the typical parameters are, but need not be, greater than 0.5mm in 

15 size in at least one dimension. Generally, in the hip joint, macro-anatomical parameters 
include thickness of the femoral shaft cortex, thickness of the femoral neck cortex, 
cortical width, hip axis length, CCD (caput-collum-diaphysis) angle, neck-shaft angle 
and width of the trochanteric region. In the spine, macro-anatomical parameters include 
thickness of the superior and inferior endplate, thickness of the anterior, lateral and 

20 posterior vertebral walls, diameter and height of the vertebral body, dimensions of the 
spinal canal and the posterior elements. 

[0054] Generally, the ray is centered onto the hip joint medial and superior to the 
greater trochanter. A calibration phantom, such as an aluminum step wedge may also 
be included in the images to calibrate gray values before further image analysis. 
25 [0055] In other embodiments, dental x-rays are preferred because of the relative 
ease and lack of expense in obtaining these images. Further, the mandible and maxilla 
are primarily composed of trabecular bone. Since the metabolic turnover of trabecular 
bone is approximately eight times greater than that of cortical bone, areas of 
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predominantly trabecular bone such as the vertebral body are preferred sites for 
measuring bone mineral density. Lang et al. (1991) Radiol Clin North Am 29:49-76. 
Thus, trabecular bone is clearly visible on the dental x-ray image, thus facilitating 
quantitative analysis of bone mineral density and structure. Jeffcoat et al. (2000) 

5 Feriodontol 23:94-1 02; Southard et al. (2000) J Dent Res 79:964-969. Further, the 
earliest bone loss in osteoporosis patients occurs in areas of trabecular bone. Multiple 
dental x-ray images are commonly made in most Americans throughout life. Indeed, 
there are approximately 750 million U.S. dental visits annually and 150 million of these 
patients result in more than 1 billion dental x-rays taken each year. Thus, the ability to 

10 diagnose osteoporosis on dental x-rays would be extremely valuable since it would 
create the opportunity for low-cost mass screening of the population. 
[0056] Preferably, x-ray imaging is performed using standard x-ray equipment, for 
instance standard dental x-ray equipment (e.g. General Electric Medical Systems, 
Milwaukee, Wl). X-rays of the incisor region and canine region are acquired using a 

15 standard x-ray imaging technique with 80 kVp and automatic exposure using a 
phototimer or using a manual technique with 10mA tube current. X-ray images are 
acquired, for example, on Kodak Ultraspeed film (Kodak, Rochester, NY). X-ray images 
may be digitized using a commercial flatbed scanner with transparency option (Acer 
ScanPremio ST). Similarly, other imaging techniques are typically performed using 

20 standard equipment, for instance, MRI or CAT equipment. 
1.1. Calibration Phantoms 
[0057] It is highly preferred that the images include accurate reference markers, 
for example calibration phantoms for assessing bone mineral density and/or bone 
structure and/or one or more macro-anatomical and/or biomechanical parameters on 

25 any given image. Calibration references (also known as calibration phantoms) for use 
in imaging technologies have been described. See, e.g., U.S. Patent No. 5,493,601 
and U.S. Patent No. 5,235,628. U.S. Patent No. 5,335,260 discloses a calibration 
phantom representative of human tissue containing variable concentrations of calcium 
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that serves as reference for quantifying calcium, bone mass and bone mineral density in 
x-ray and CT imaging systems. However, currently available calibration phantoms are 
not always accurate. Because bone mineral density accounts for considerably less than 
100% of fracture risk in osteoporosis (Ouyang et al. (1997) Calif Tissue Int, 60:139-147) 

5 some of the methods and devices described herein are designed to assess not only 
bone mineral density but also bone structure and, in addition, macro-anatomical and/or 
biomechanical parameters. By assessing two or more of these parameters, more 
accurate testing and screening can be provided for conditions such as osteoporosis. 
[0058] Thus, in certain aspects, the current invention provides for methods and 

10 devices that allow accurate quantitative assessment of information contained in an x-ray 
such as density of an anatomic structure and/or morphology of an anatomic structure. 
Any suitable calibration phantom can be used, for example, one that comprises 
aluminum or other radio-opaque materials. U.S. Patent No. 5,335,260 describes other 
calibration phantoms suitable for use in assessing bone mineral density in images. 

15 Examples of other suitable calibration reference materials can be fluid or fluid-like 
materials, for example, one or more chambers filled with varying concentrations of 
calcium chloride or the like. 

[0059] Numerous calibration phantoms (or reference calibrations) can be used in 
the practice of the present invention. Typically, the system used to monitor bone 

20 mineral density and/or bone structure and/or one or more macro-anatomical and/or 
biomechanical parameters in a target organism comprises an image (e.g., a dental or 
hip radiograph), which provides information on the subject; an assembly including a 
calibration phantom, which acts as a reference for the data in the image; and at least 
one data processing system, which evaluates and processes the data from the image 

25 and/or from the calibration phantom assembly. 

[0060] It will be readily apparent that a calibration phantom can contain a single, 
known density or structure reference. Furthermore, a gradient in density can be 
achieved by varying the thickness or the geometry of the calibration phantom along the 
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path of the x-ray beam, for example, by using a V-shape of the calibration phantom of 
varying thickness (Fig. 4). The calibration phantom can also include angles. For 
example, the calibration phantom can be "T"-shaped or "L"-shaped thereby including 
one or more 90 degree angles. 
5 [00@1] The calibration phantom can contain several different areas of different 
radio-opacity. For example, the calibration phantom can have a step-like design, 
whereby changes in local thickness of the wedge result in differences in radio-opacity. 
Stepwedges using material of varying thickness are frequently used in radiology for 
quality control testing of x-ray beam properties. By varying the thickness of the steps, 
10 the intensity and spectral content of the x-ray beam in the projection image can be 
varied. Stepwedges are commonly made of aluminum, copper and other convenient 
and homogeneous materials of known x-ray attenuation properties. Stepwedge-like 
phantoms can also contain calcium phosphate powder or calcium phosphate powder in 
molten paraffin. 

15 [0062] Alternatively, continuous wedges may be used or the calibration reference 
may be designed such that the change in radio-opacity is from periphery to center (for 
example in a round, ellipsoid, rectangular, triangular of other shaped structure). As 
noted above, the calibration reference can also be constructed as plurality of separate 
chambers, for example fluid filled chambers, each including a specific concentration of a 

20 reference fluid (e.g., calcium chloride). In addition to one or more fluids, a calibration 
phantom can also contain metal powder, e.g. aluminum or steel powder, embedded 
within it (for example, embedded in a plastic). 

[0063] In certain embodiments, the calibration phantom is specifically designed to 
serve as a reference for bone structure (e.g., trabecular spacing, thickness and the like). 
25 For example, the calibration wedge can contain one or more geometric patterns with 
known dimensions, e.g. a grid whereby the spacing of a grid, thickness of individual grid 
elements, etc. are known. This known geometric pattern of radio-opaque elements in 
the calibration phantom can be used to improve the accuracy of measurements of 
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trabecular bone structure in an x-ray. Such measurements of trabecular bone structure 
can include, but are not limited to, trabecular spacing, trabecular length and trabecular 
thickness. Such measurements of trabecular spacing, trabecular length and trabecular 

thickness can, for example, be performed in a denial or spine or hip x-ray. These 
5 calibration phantoms can be made up of a variety of materials include, plastics, metals 
and combinations thereof. Further, the reference components can be solid, powdered, 
fluid or combinations thereof. Thus, the calibration wedge can also be used to improve 
measurements of bone structure. 

[0064] In certain embodiments, the calibration phantom is specifically designed to 

10 serve as a reference for macro-anatomical parameters (e.g., in the hip joint, thickness of 
the femoral shaft cortex, thickness of the femoral neck cortex, cortical width, hip axis 
length, CCD (caput-collum-diaphysis) angle, neck-shaft angle and width of the 
trochanteric region; and in the spine, thickness of the superior and inferior endplate, 
thickness of the anterior, lateral and posterior vertebral walls, diameter and height of the 

15 vertebral body, dimensions of the spinal canal and the posterior elements). For 

example, the calibration wedge can contain one or more geometric patterns with known 
dimensions, e.g. a grid whereby the spacing of a grid, thickness of individual grid 
elements, etc. are known. This known geometric pattern of radio-opaque elements in 
the calibration phantom can be used to improve the accuracy of measurements of 

20 macro-anatomical and/or biomechanical parameters in an x-ray, for example by aiding 
in the correction of image magnification. Such measurements of macro-anatomical 
parameters can, for example, be performed in a dental or spine or hip x-ray. These 
calibration phantoms can be made up of a variety of materials include, plastics, metals 
and combinations thereof. Further, the reference components can be solid, powdered, 

25 fluid or combinations thereof. Thus, the calibration wedge can also be used to improve 
measurements of bone structure. 

[0065] Since the present invention contemplates analysis of dental x-ray images 
for information on bone structure, bone mineral density or both structure and density, it 
will be apparent that calibration phantoms will be selected based on whether structure, 
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density or both are being measured. Thus, one or more calibration phantoms may be 
present. 

[0068] Whatever the overall shape or composition of the calibration phantom, 
when present, the at least one marker be positioned at a known density and/or structure 

5 in the phantom. Furthermore, it is preferred that at least one geometric shape or pattern 
is included in the calibration phantom. Any shape can be used including, but not limited 
to, squares, circles, ovals, rectangles, stars, crescents, multiple-sided objects (e.g., 
octagons), V- or U-shaped, inverted V- or U-shaped, irregular shapes or the like, so 
long as their position is known to correlate with a particular density of the calibration 

10 phantom. In preferred embodiments, the calibration phantoms described herein are 
used in 2D planar x-ray imaging. 

[0067] The calibration phantoms can be imaged before or after the x-ray image is 
taken. Alternatively, the calibration phantom can be imaged at the same time as the x- 
ray image. The calibration phantom can be physically connected to an x-ray film and/or 

15 film holder. Such physical connection can be achieved using any suitable mechanical 
or other attachment mechanism, including but not limited to adhesive, a chemical bond, 
use of screws or nails, welding, a Velcro™ strap or Velcro™ material and the like. 
Similarly, a calibration phantom can be physically connected to a detector system or a 
storage plate for digital x-ray imaging using one or more attachment mechanisms (e.g., 

20 a mechanical connection device, a Velcro™ strap or other Velcro™ material, a chemical 
bond, use of screws or nails, welding and an adhesive). The external standard and the 
film can be connected with use of a holding device, for example using press fit for both 
film and external standard. 

[0068] Additionally, the calibration phantom assembly can be attached to an 
25 anatomical structure, for example one or more teeth, mucus membranes, the mandible 
and/or maxilla. For instance, the calibration phantom can be attached (e.g., via 
adhesive attachment means) to the epithelium or mucous membrane inside overlying 
the mandible or the maxilla. Alternatively, the calibration phantom can be placed on or 
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adjacent to a tooth, for example, a V- or U-shaped (in the case of the maxilla) or an 
inverted V- or U-shaped (in the case of the mandible) calibration phantom can be used. 
The opening of the V or U will be in contact with the free edge of at least one tooth or 
possibly several teeth (Fig. 4). 

5 [0069] In preferred embodiments, when an x-ray of an anatomic structure or a 
non-living object is acquired a calibration phantom is included in the field of view. Any 
suitable calibration phantom can be used, for example, one that comprises aluminum or 
other radio-opaque materials. U.S. Patent No. 5,335,260 describes other calibration 
phantoms suitable for use in assessing bone mineral density in images. Examples of 

10 other suitable calibration reference materials can be fluid or fluid-like materials, for 

example, one or more chambers filled with varying concentrations of calcium chloride or 
the like. In a preferred embodiment, the material of the phantom is stainless steel {e.g., 
AISI grade 316 comprising carbon (0.08%); manganese (2%); silicon (1%); phosphorus 
(0.045%); sulphur (0.03%); nickel (10-14%); chromium (16-18%); molybdenum (2-3%); 

15 plus iron to make up 100%). The relative percentages of the components may be with 
respect to weight or volume. 

[0070] It will be apparent that calibration phantoms suitable for attachment to an 
anatomical structure can have different shapes depending on the shape of the 
anatomical structure (e.g., tooth or teeth) on which or adjacent to which it will be placed 

20 including, but not limited to, U-shaped, V-shaped, curved, flat or combinations thereof. 
For example, U-shaped (or inverted U-shaped) calibration phantoms can be positioned 
on top of molars while V-shaped (or inverted V-shaped) calibration phantoms can be 
positioned on top of incisors. Further, it will be apparent that in certain instances (e.g., 
teeth on the mandible), the calibration phantom can rest on top of the tooth just based 

25 on its gravity or it can be attached to the tooth (e.g., using adhesive). In the case of the 
teeth on the maxilla, the calibration phantom will typically be attached to the tooth, for 
example with use of an adhesive. 
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[0071] Any of these attachments may be permanent or temporary and the 
calibration phantom can be integral (e.g., built-in) to the film, film holder and/or detector 
system or can be attached or positioned permanently or temporarily appropriately after 
the film and/or film holder is produced. Thus, the calibration phantom can be designed 

5 for single-use (e.g., disposable) or for multiple uses with different x-ray images. Thus, 
in certain embodiments, the calibration phantom is reusable and, additionally, can be 
sterilized between uses. Integration of a calibration phantom can be achieved by 
including a material of known x-ray density between two of the physical layers of the x- 
ray film. Integration can also be achieved by including a material of known x-ray density 

10 within one of the physical layers of the x-ray film. Additionally, the calibration phantom 
can be integrated into the film cover. A calibration phantom or an external standard can 
also be integrated into a detector system or a storage plate for digital x-ray imaging. For 
example, integration can be achieved by including a material of known x-ray density 
between two of the physical layers of the detector system or the storage plate. 

15 Integration can also be achieved by including a material of know x-ray density within 
one of the physical layers of the detector system or the storage plate. 
[0072] In certain embodiments, for example those embodiments in which the 
calibration phantom is temporarily attached to a component of the x-ray assembly 
system (e.g., x-ray film holder, x-ray film, detector system or the like), cross-hairs, lines 

20 or other markers may be placed on the apparatus as indicators for positioning of the 
calibration phantom. These indicators can help to ensure that the calibration phantom 
is positioned such that it doesn't project on materials that will alter the apparent density 
in the resulting image. 

[0073] Any of the calibration phantom-containing assemblies described herein 
25 can be used in methods of analyzing and/or quantifying bone structure and/or one or 
more macro-anatomical and/or biomechanical parameters (or bone mineral density) in 
an x-ray image. The methods generally involve simultaneously imaging or scanning the 
calibration phantom and another material (e.g., bone tissue from a subject) for the 
purpose of quantifying the density of the imaged material (e.g., bone mass). In the case 
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of dental radiographs, the calibration phantom, the x-ray tube or dental x-ray film is 
typically positioned in a manner to ensure inclusion of the calibration phantom and a 
portion of the mandible and/or maxilla on the dental x-ray image. Preferably, the 
calibration phantom, the x-ray tube and the dental x-ray film are positioned so that at 

5 least a portion of the section of the mandible or maxilla included on the image will 
contain predominantly trabecular bone rather than cortical bone. 
[00741 Thus, under the method of the present invention, the calibration phantom 
is preferably imaged or scanned simultaneously with the individual subject, although the 
invention allows for non-simultaneous scanning of the phantom and the subject. 

10 Methods of scanning and imaging structures by x-ray imaging technique are well 

known. By placing the calibration phantom in the x-ray beam with the subject, reference 
calibration samples allow corrections and calibration of the absorption properties of 
bone. When the phantom is imaged or scanned simultaneously with each subject, the 
variation in x-ray beam energy and beam hardening are corrected since the phantom 

15 and the subject both see the same x-ray beam spectrum. Each subject, having a 
different size, thickness, muscle-to-fat ratio, and bone content, attenuate the beam 
differently and thus change the effective x-ray beam spectrum. It is necessary that the 
bone-equivalent calibration phantom be present in the same beam spectrum as the 
subject's bone to allow accurate calibration. 

20 [0075] X-ray imaging assemblies that are currently in use do not take into 
account the position of the calibration phantom in relation to the structures being 
imaged. Thus, when included in known assemblies, calibration phantom(s) are often 
positioned such that they project on materials or structures {e.g., bone) that alter 
apparent density of the calibration phantom in the resulting x-ray image. Clearly, this 

25 alteration in apparent density will affect the accuracy of the calibration phantom as a 
reference for determining bone mineral density, structure or macro-anatomical 
parameters. Therefore, it is an object of the invention to provide methods in which the 
calibration phantom projects free of materials or structures that will alter the apparent 
density of the reference. In the context of dental x-rays, for instance, the methods 
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described herein ensure that the calibration phantom projects free of bone (e.g., teeth, 
jaw) tissue. This can be accomplished in a variety of ways, for example, positioning the 
calibration phantom in the x-ray film or in the x-ray film holder such that it will appear 
between the teeth in the dental x-ray. 

5 [D07@] The calibration phantom materials and methods of the present invention 
are preferably configured to be small enough and thin enough to be placed inside the 
mouth, and the method of the present invention can be used to quantify bone mass 
using standard dental x-ray systems, for example by including temporary or permanent 
calibration phantoms in dental x-ray film holders. Further, it is highly desirable that the 

10 calibration phantom be positioned so that at least a portion doesn't project on structures 
or materials that will alter the apparent density or structural characteristics of the 
calibration phantoms. It is also preferable to position calibration phantom at a defined 
distance relative to at least one tooth or the mandible or the maxilla whereby a 
substantial portion of the calibration phantom projects free of said tooth, said mandible 

15 or said maxilla on the x-ray image. Any suitable distance can be used, for example 
between about 1 mm and 5 cm or any value therebetween. 

[0077] A cross-calibration phantom can be used to optimize system performance, 
e.g. x-ray tube settings or film processor settings, or to improve the comparability of 
different machines or systems, typically located at different sites. For this purpose, a 

20 separate image may be obtained which does not include a patient or a body part. The 
image includes the primary calibration phantom used in patients, e.g. a step-wedge of 
known density, and the cross-calibration phantom. The apparent density of the primary 
calibration phantom is then calibrated against the density of the cross-calibration 
phantom. The resultant cross-calibration of the primary phantom can help to improve 

25 the accuracy of measurements of bone density, bone structure and macro-anatomical 
and/or biomechanical parameters. It can also help improve the overall reproducibility of 
the measurements. In one embodiment of the invention, an x-ray technologist or a 
dental hygienist will perform a cross-calibration test once a day, typically early in the 
morning, prior to the first patient scans. The results of the cross-calibration or the entire 
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cross-calibration study can be transmitted via a network to a central computer. The 
central computer can then perform adjustments designed to maintain a high level of 

comparability between different systems. 

1 .2. Inherent Reference barkers 

5 [0078] In certain embodiments of the invention, information inherent in the 
anatomic structure or the non-living object can be used to estimate the density and/or 
structure and/or macro-anatomy of selected bone regions of interest within the anatomic 
structure or the non-living object. For example, since the density of muscle, fat, water 
(e.g., soft tissue), metal (e.g., dental fillings) and air are typically known, the density of 

10 air surrounding an anatomic structure or non-living object, the density of subcutaneous 
fat, and the density of muscle tissue can be used to estimate the density of a selected 
region of bone, for example within the distal radius. For instance, a weighted mean can 
be determined between one or more of the internal standards (e.g., air, water, metal, 
and/or fat) and used as internal standards to determine bone density in the same x-ray 

15 image. Similarly, the density of a tooth or a portion of a tooth can be used to estimate 
the density of a selected region of bone, e.g. an area in the mandible. 
[0079] The information inherent in said anatomic structure can also be combined 
with information provided by the calibration phantom and the combination can result in 
an improved accuracy of the calibration phantom. 

20 1.3. Holders and Hygienic Covers 

[0080] As noted above, in certain embodiments, a holder can be used to position 
the calibration phantom. The holder can be U-shaped or V-shaped (Fig. 5) for ease in 
attachment to a tooth. The attachment can be, for example, with an adhesive. The 
calibration phantom, in turn, can be attached to the holder. Similarly, the calibration 

25 phantom can be attached to holders comprising one or more molds of at least one or 
more teeth. Additionally, the holder can be used to position both the film and the 
calibration phantom relative to the osseous structure that will be included in the x-ray 
image. In another embodiment, a holding device that can hold the x-ray film is 
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integrated in the calibration phantom. This holding device can hold the film in place 
prior to taking the x-ray. The holding device can be spring-loaded or use other means 
such as mechanical means of holding and stabilizing the x-ray film. 
[0081] In certain embodiments, the holder may comprise a disposable or 

5 sterilizeable hygienic cover. See, e.g., WO 99/08598, the disclosure of which is 

incorporated by reference herein in its entirety. Furthermore, the holder may comprise 
multiple components, for example, the calibration phantom and a integrated or 
insertable bolus back that can serve to enhance the accuracy of the calibration phantom 
by accounting for the effect of soft tissue that may project with the calibration phantom 

10 and/or with the bone. 

[0082] In certain embodiments, the calibration phantom can be configured so that 
it stabilizes against the surrounding tissues on its own without the use of an additional 
holder, The calibration phantom can be protected with a hygienic cover. 
[0083] The holder (e.g., hygienic cover) may be comprised of a rigid material, a 

15 flexible material or combinations thereof. Furthermore, the holder may include one or 
more pockets/compartments adapted to receive additional components such as the 
calibration phantom, a bolus back or the like. Additionally, one or more portions of the 
holder may be radiolucent. 

2.0. Analysis and Manipulation of Data 

20 [0084] The data obtained from images taken as described above is then 

preferably analyzed and manipulated. Thus, the systems and assemblies described 
herein can also include one or more computational units designed, for example, to 
analyze bone density or bone structure or macro-anatomical and/or biomechanical data 
in the image; to identify an anatomical landmark in an anatomical region; to correct for 

25 soft tissue measurements; and/or to evaluate bone density and structure and macro- 
anatomy of the image. As will be appreciated by those of skill in the art, the 
computational unit can include any software, chip or other device used for calculations. 
Additionally, the computational unit may be designed to control the imaging assembly or 
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detector (as well as other parameters related to the detector(s)). Other applications of 
the computational unit to the methods and devices described herein will be recognized 
by those skilled in the art. The computational unit may be used for any other application 
related to this technology that may be facilitated with uss of computer software or 

5 hardware. The computational unit can also further comprise a database comprising, for 
example, reference anatomical maps and the computational unit is further designed to 
compare the anatomical map with the reference anatomical map. The reference 
anatomical map may be historic (from the same or another patient, generated as part of 
an interrogation protocol), or theoretical or any other type of desired reference map. 

10 [0085] Any image can be analyzed in order to obtain and manipulate data. Thus, 
data points, derived data, and data attributes database according to the present 
invention may comprise the following: (1) the collection of data points, said data points 
comprising information obtained from an image, for example, bone mineral density 
information or information on bone structure (architecture); and (2) the association of 

15 those data points with relevant data point attributes. The method may further comprise 
(3) determining derived data points from one or more direct data points and (4) 
associating those data points with relevant data point attributes. The method may also 
comprise (5) collection of data points using a remote computer whereby said remote 
computer operates in a network environment. 

20 [0086] In certain preferred embodiments, the information is obtained from a 

dental x-ray image. As described herein, dental x-ray images can be acquired at a local 
site using known techniques. If the x-ray image was captured using conventional x-ray 
film, the data points (information) of the x-ray image can be digitized using a scanning 
device. The digitized x-ray image information can then be transmitted over the network, 

25 e.g. the Internet, into a remote computer or server. If the x-ray image was acquired 

using digital acquisition techniques, e.g. using phosphorus plate systems or selenium or 
silicon detector systems, the x-ray image information is already available in digital 
format. In this case the image can be transmitted directly over the network, e.g. the 
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Internet. The information can also be compressed and/or encrypted prior to 
transmission. Transmission can also be by other methods such as fax, mail or the like. 

2.1. Data Points 

[0087] Thus, the methods of and compositions described herein make use of 

5 collections of data sets of measurement values, for example measurements of bone 
structure and/or bone mineral density from x-ray images. Records may be formulated in 
spreadsheet-like format, for example including data attributes such as date of x-ray, 
patient age, sex, weight, current medications, geographic location, etc. The database 
formulations may further comprise the calculation of derived or calculated data points 

10 from one or more acquired data points. A variety of derived data points may be useful 
in providing information about individuals or groups during subsequent database 
manipulation, and are therefore typically included during database formulation. Derived 
data points include, but are not limited to the following: (1) maximum bone mineral 
density, determined for a selected region of bone or in multiple samples from the same 

15 or different subjects; (2) minimum bone mineral density, determined for a selected 
region of bone or in multiple samples from the same or different subjects; (3) mean 
bone mineral density, determined for a selected region of bone or in multiple samples 
from the same or different subjects; (4) the number of measurements that are 
abnormally high or low, determined by comparing a given measurement data point with 

20 a selected value; and the like. Other derived data points include, but are not limited to 
the following: (1) maximum value of a selected bone structure parameter, determined 
for a selected region of bone or in multiple samples from the same or different subjects; 
(2) minimum value of a selected bone structure parameter, determined for a selected 
region of bone or in multiple samples from the same or different subjects; (3) mean 

25 value of a selected bone structure parameter, determined for a selected region of bone 
or in multiple samples from the same or different subjects; (4) the number of bone 
structure measurements that are abnormally high or low, determined by comparing a 
given measurement data point with a selected value; and the like. Other derived data 
points include, but are not limited to the following: (1) maximum value of a selected 



25 



WO 2004/086972 



PCT7US2004/00<M<>5 



macro-anatomical and/or biomechanical parameter, determined for a selected region of 
bone or in multiple samples from the same or different subjects; (2) minimum value of a 
selected macro-anatomical and/or biomechanical parameter, determined for a selected 
region of bone or in multiple samples from the same or differeni subjects; (3) mean 
5 value of a selected macro-anatomical and/or biomechanical parameter, determined for 
a selected region of bone or in multiple samples from the same or different subjects; (4) 
the number of macro-anatomical and/or biomechanical measurements that are 
abnormally high or low, determined by comparing a given measurement data point with 
a selected value; and the like. Other derived data points will be apparent to persons of 

10 ordinary skill in the art in light of the teachings of the present specification. The amount 
of available data and data derived from (or arrived at through analysis of) the original 
data provide provides an unprecedented amount of information that is very relevant to 
management of bone related diseases such as osteoporosis. For example, by 
examining subjects over time, the efficacy of medications can be assessed. 

15 [0088] Measurements and derived data points are collected and calculated, 
respectively, and may be associated with one or more data attributes to form a 
database. The amount of available data and data derived from (or arrived at through 
analysis of) the original data provide provides an unprecedented amount of information 
that is very relevant to management of bone related diseases such as osteoporosis. 

20 For example, by examining subjects over time, the efficacy of medications can be 
assessed. 

[0089] Data attributes can be automatically input with the x-ray image and can 
include, for example, chronological information (e.g., DATE and TIME). Other such 
attributes may include, but are not limited to, the type of x-ray imager used, scanning 
25 information, digitizing information and the like. Alternatively, data attributes can be input 
by the subject and/or operator, for example subject identifiers, i.e. characteristics 
associated with a particular subject. These identifiers include but are not limited to the 
following: (1) a subject code (e.g., a numeric or alpha-numeric sequence); (2) 
demographic information such as race, gender and age; (3) physical characteristics 
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such as weight, height and body mass index (BMI); (4) selected aspects of the subject's 
medical history (e.g., disease states or conditions, etc.); and (5) disease-associated 
characteristics such as the type of bone disorder, if any; the type of medication used by 
the subject. In the practice of the present invention, each data point would typically be 
5 identified with the particular subject, as well as the demographic, etc. characteristic of 
that subject. 

[0090] Other data attributes will be apparent to persons of ordinary skill in the art 
in light of the teachings of the present specification. 

2.2. Storage of Data Sets and Association of Data Points with 
10 Relevant Data Attributes 

[0091] A number of formats exist for storing data sets and simultaneously 
associating related attributes, including but not limited to (1) tabular, (2) relational, and 
(3) dimensional. In general the databases comprise data points, a numeric value which 
correspond to physical measurement (an "acquired" datum or data point) or to a single 
15 numeric result calculated or derived from one or more acquired data points that are 
obtained using the various methods disclosed herein. The databases can include raw 
data or can also include additional related information, for example data tags also 
referred to as "attributes" of a data point. The databases can take a number of different 
forms or be structured in a variety of ways. 

20 [0092] The most familiar format is tabular, commonly referred to as a 

spreadsheet. A variety of spreadsheet programs are currently in existence, and are 
typically employed in the practice of the present invention, including but not limited to 
Microsoft Excel spreadsheet software and Corel Quattro spreadsheet software. In this 
format, association of data points with related attributes occurs by entering a data point 

25 and attributes related to that data point in a unique row at the time the measurement 
occurs. 

[0093] Further, rational, relational (Database Design for Mere Mortals, by Michael 
J. Hernandez, 1997, Addison-Wesley Pub. Co., publisher; Database Design for 
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Smarties, by Robert J. Muller, 1999, Morgan Kaufmann Publishers, publisher; 
Relational Database Design Clearly Explained, by Jan L. Harrington, 1998, Morgan 
Kaufmann Publishers, publisher) and dimensional (Data-Parallel Computing, by V.B. 
IWiuchnick, etal., 1996, International Thomson Publishing, publisher; Understanding 
5 Fourth Dimensions, by David Graves, 1993, Computerized Pricing Systems, publisher) 
database systems and management may be employed as well. 
[0094] Relational databases typically support a set of operations defined by 
relational algebra. Such databases typically include tables composed of columns and 
rows for the data included in the database. Each table of the database has a primary 

10 key, which can be any column or set of columns, the values for which uniquely identify 
the rows in a table. The tables in the database can also include a foreign key that is a 
column or set of columns, the values of which match the primary key values of another 
table. Typically, relational databases also support a set of operations (e.g., select, join 
and combine) that form the basis of the relational algebra governing relations within the 

15 database. 

[0095] Such relational databases can be implemented in various ways. For 
instance, in Sybase® (Sybase Systems, Emeryville, CA) databases, the tables can be 
physically segregated into different databases. With Oracle® (Oracle Inc., Redwood 
Shores, CA) databases, in contrast, the various tables are not physically separated, 
20 because there is one instance of work space with different ownership specified for 

different tables. In some configurations, databases are all located in a single database 
(e.g., a data warehouse) on a single computer. In other instances, various databases 
are split between different computers. 

[0096] It should be understood, of course, that the databases are not limited to 
25 the foregoing arrangements or structures. A variety of other arrangements will be 
apparent to those of skill in the art. 

2.3. Data Manipulation 
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[0097] Data obtained from x-ray images as described herein can be manipulated, 
for example, using a variety of statistical analyses, to produce useful information. The 
databases of the present invention may be generated, for example, from data collected 
for an individual or from a selected group of individuals over a defined period of time 
5 (e.g., days, months or years), from derived data, and from data attributes. 

[0098] For example, data may be aggregated, sorted, selected, sifted, clustered 
and segregated by means of the attributes associated with the data points. A number of 
data mining software programs exist which may be used to perform the desired 
manipulations. 

10 [0099] Relationships in various data can be directly queried and/or the data 

analyzed by statistical methods to evaluate the information obtained from manipulating 
the database. 

[0100] For example, a distribution curve can be established for a selected data 
set, and the mean, median and mode calculated therefor. Further, data spread 

15 characteristics, e.g. variability, quartiles and standard deviations can be calculated. 
[0101] The nature of the relationship between any variables of interest can be 
examined by calculating correlation coefficients. Useful methods for doing so include 
but are not limited to the following: Pearson Product Moment Correlation and Spearman 
Rank Order Correlation. 

20 [0102] Analysis of variance permits testing of differences among sample groups 
to determine whether a selected variable has a discernible effect on the parameter 
being measured. 

[0103] Non-parametric tests may be used as a means of testing whether 
variations between empirical data and experimental expectancies are attributable 
25 merely to chance or to the variable or variables being examined. These include the Chi 
Square test, the Chi Square Goodness of Fit, the 2 x 2 Contingency Table, the Sign 
Test, and the Phi Correlation Coefficient. 
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[0104] There are numerous tools and analyses available in standard data mining 
software that can be applied to the analysis of the databases of the present invention. 
Such tools and analyses include, but are not limited to, cluster analysis, factor analysis, 
decision trees, neural networks, rule induction, data driven modeling, and data 
5 visualization. Some of the more complex methods of data mining techniques are used 
to discover relationships that are more empirical and data-driven, as opposed to theory- 
driven, relationships. 

[01 05] Exemplary data mining software that can be used in analysis and/or 
generation of the databases of the present invention includes, but is not limited to: Link 

10 Analysis (e.g., Associations analysis, Sequential Patterns, Sequential time patterns and 
Bayes Networks); Classification (e.g., Neural Networks Classification, Bayesian 
Classification, k-nearest neighbors classification, linear discriminant analysis, Memory 
based Reasoning, and Classification by Associations); Clustering (e.g., k-Means 
Clustering, demographic clustering, relational analysis, and Neural Networks 

15 Clustering); Statistical methods (e.g., Means, Std dev, Frequencies, Linear Regression, 
non-linear regression, t-tests, F-test, Chi2 tests, Principal Component Analysis, and 
Factor Analysis); Prediction (e.g., Neural Networks Prediction Models, Radial Based 
Functions predictions, Fuzzy logic predictions, Times Series Analysis, and Memory 
based Reasoning); Operating Systems; and Others (e.g., Parallel Scalability, Simple 

20 Query Language functions, and C++ objects generated for applications). Companies 
that provide such software include, for example, the following: Adaptative Methods 
Group at UTS (UTS City Campus, Sydney, NSW 2000), CSI®, Inc., (Computer Science 
Innovations, Inc. Melbourne, Florida), IBM® (International Business Machines 
Corporation, Armonk, NY), Oracle® (Oracle Inc., Redwood Shores, CA) and SAS® 

25 (SAS Institute Inc., Cary, NC). 

[0106] These methods and processes may be applied to the data obtained using 
the methods described herein, for example, databases comprising, x-ray image data 
sets, derived data, and data attributes. 
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[0107] In certain embodiments, data (e.g., bone structural information or macro- 
anatomical and/or biomechanical information or bone mineral density information) is 
obtained from normal control subjects using the methods described herein. These 
databases are typically referred to as "reference databases" and can be used to aid 

5 analysis of any given subject's x-ray image, for example, by comparing the information 
obtained from the subject to the reference database. Generally, the information 
obtained from the normal control subjects will be averaged or otherwise statistically 
manipulated to provide a range of "normal" (reference) measurements. Suitable 
statistical manipulations and/or evaluations will be apparent to those of skill in the art in 

10 view of the teachings herein. 

[0108] ADD z-score, T-scores here; other statistical measurements that you 
deem important. The comparison of the subject's x-ray information to the reference 
database can be used to determine if the subject's bone information falls outside the 
normal range found in the reference database or is statistically significantly different 

15 from a normal control. Data comparison and statistical significance can be readily 

determined by those of skill in the art using for example the z-test or t-test statistics for 
continuous variables, the chi-square test or Fisher's exact test for categorical data and 
the rank-sum test or Kruskal-Wallis test for ranked data. The use of reference 
databases in the analysis of x-ray images facilitates that diagnosis, treatment and 

20 monitoring of bone conditions such as osteoporosis. 

[0109] For a general discussion of statistical methods applied to data analysis, 
see Applied Statistics for Science and Industry, by A. Romano, 1977, Allyn and Bacon, 
publisher. 

[01 1 0] The data is preferably stored and manipulated using one or more 
25 computer programs or computer systems. These systems will typically have data 
storage capability (e.g., disk drives, tape storage, CD-ROMs, etc.). Further, the 
computer systems may be networked or may be stand-alone systems. If networked, the 
computer system would be able to transfer data to any device connected to the 
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networked computer system for example a medical doctor or medical care facility using 
standard e-mail software, a central database using database query and update software 
(e.g., a data warehouse of data points, derived data, and data attributes obtained from a 
large number of subjects). Alternatively, a user could access from a doctor's office or 
5 medical facility, using any computer system with Internet access, to review historical 
data that may be useful for determining treatment. 

[0111] If the networked computer system includes a World Wide Web application, 
the application includes the executable code required to generate database language 
statements, for example, SQL statements. Such executables typically include 

10 embedded SQL statements. The application further includes a configuration file that 
contains pointers and addresses to the various software entities that are located on the 
database server in addition to the different external and internal databases that are 
accessed in response to a user request. The configuration file also directs requests for 
database server resources to the appropriate hardware, as may be necessary if the 

15 database server is distributed over two or more different computers. 

[01 1 2] Usually each networked computer system includes a World Wide Web 
browser that provides a user interface to the networked database server. The 
networked computer system is able to construct search requests for retrieving 
information from a database via a Web browser. With access to a Web browser users 

20 can typically point and click to user interface elements such as buttons, pull down 

menus, and other graphical user interface elements to prepare and submit a query that 
extracts the relevant information from the database. Requests formulated in this 
manner are subsequently transmitted to the Web application that formats the requests 
to produce a query that can be used to extract the relevant information from the 

25 database. 

[0113] When Web-based applications are utilized, the Web application accesses 
data from a database by constructing a query in a database language such as Sybase 
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or Oracle SQL which is then transferred to a relational database management system 
that in turn processes the query to obtain the pertinent information from the database. 
[0114] Accordingly, in one aspect the present invention describes* a method of 
providing data obtained from x-ray images on a network, for example 'the Internet, and 

5 methods of using this connection to provide real-time and delayed data analysis. The 
central network can also allow access by the physician to a subject's data. Similarly, an 
alert could be sent to the physician if a subject's readings are out of a predetermined 
range, etc. The physician can then send advice back to the patient via e-mail or a 
message on a web page interface. Further, access to the entire database of data from 

10 all subjects may be useful for statistical or research purposes. Appropriate network 
security features (e.g., for data transfer, inquiries, device updates, etc.) are of course 
employed. 

[0115] Further, a remote computer can be used to analyze the x-ray that has 
been transmitted over the network automatically. For example, x-ray density information 
15 or structural information about an object can be generated in this fashion. X-ray density 
information can, for example, be bone mineral density. If used in this fashion, the test 
can be used to diagnose bone-related conditions such as osteoporosis. 

2.4. Graphical User Interface 

[01 16] In certain of the computer systems, an interface such as an interface 
20 screen that includes a suite of functions is included to enable users to easily access the 
information they seek from the methods and databases of the invention. Such 
interfaces usually include a main menu page from which a user can initiate a variety of 
different types of analyses. For example, the main menu page for the databases 
generally include buttons for accessing certain types of information, including, but not 
25 limited to, project information, inter-project comparisons, times of day, events, dates, 
times, ranges of values, etc. 

2.5. Computer Program Products 
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[01 17] A variety of computer program products can be utilized for conducting the 
various methods and analyses disclosed herein. In general, the computer program 
products comprise a computer-readable medium and the code necessary to perform the 
methods set forth supra. The computer-readable medium on which the program 
5 instructions are encoded can be any of a variety of known medium types, including, but 
not limited to, microprocessors, floppy disks, hard drives, ZIP drives, WORM drives, 
magnetic tape and optical medium such as CD-ROMs. 

[Q I i 8] For example, once an image or data from that image is transmitted via a local 
or long-distance computer network and the data received by a remote computer or a 

10 computer connected to the remote network computer, an analysis of the morphology 
and density of the bone can be performed, for example using suitable computer 
programs. This analysis of the object's morphology can occur in two-dimensions or 
three-dimensions. For example, in imaging osseous structures, such analysis of the 
transmitted x-ray image can be used to measure parameters that are indicative or 

15 suggestive of bone loss or metabolic bone disease. Such parameters include all current 
and future parameters that can be used to evaluate osseous structures. For example, 
such parameters include, but are not limited to, trabecular spacing, trabecular thickness, 
trabecular connectivity and intertrabecular space. 

[0119] Information on the morphology or 2D or 3D structure of an anatomic object 
20 can be derived more accurately, when image acquisition parameters such as spatial 
resolution are known. Other parameters such as the degree of cone beam distortion can 
also be helpful in this setting. 

[0120] As noted above, an image can be transmitted from a local site into a 
remote server and the remote server can perform an automated analysis of the image. 
25 Further, the remote server or a computer connected to the remote server can then 

generate a diagnostic report. Thus, in certain embodiments, a computer program (e.g., 
on the remote server or on a computer connected to the remote server) can generate 
charges for the diagnostic report. The remote server can then transmit the diagnostic 
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report to a physician, typically the physician who ordered the test or who manages the 
patient. The diagnostic report can also be transmitted to third parties, e.g. health 
insurance companies. Such transmission of the diagnostic report can occur 
electronically (e.g. via e-mail), via mail, fax or other means of communication. All or 
5 some of the transmitted information (e.g., patient identifying information) can be 
encrypted to preserve confidentiality of medical records. 

[0121] Thus, one exemplary system is described herein for analyzing bone 
morphology or structure in a subject system via a dental x-ray that includes at least a 
portion of the mandible and/or maxilla of a subject, followed by evaluation or the x-ray 
10 image. Dental x-rays are obtained in any conventional method. The x-ray produces an 
image that can be interpreted (for example, employing a selected algorithm and/or 
computer program) by an associated system controller to provide a bone mineral 
density or bone structure evaluation for display. 

[0122] In a further aspect of the present invention, the monitoring system can 
15 comprise two or more components, in which a first component comprises an x-ray 

image and calibration phantom that are used to extract and detect bone-related data on 
the subject, and a second component that receives the data from the first component, 
conducts data processing on the data and then displays the processed data. 
Microprocessor functions can be found in one or both components. The second 
20 component of the monitoring system can assume many forms 
3.0.0.0 Correction Factors 

[0123] Although the presence of calibration phantoms greatly aids in increasing 
the accuracy of data obtained from images such as dental, hip or spine x-rays, the 
present inventors also recognize that, in certain instances, there may be a need to apply 
25 one or more correction factors to further enhance accuracy of the data obtained from 
any given x-ray image. Such correction factors will take into account one or more of a 
wide variety of influences (e.g.,. soft tissue thickness, region from which the data is 
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extracted and the like) that can alter apparent density or structure information on the 
image. 

[01 24] In this regard, one or more reference databases can be used for 
calibration and normaliiation purposes. For example, image normalization or correction 
5 of soft-tissue attenuation can be performed using patient characteristic data such as 
patient weight, height and body mass index. In one example, a higher soft-tissue 
attenuation can be assumed in high weight and low height subjects; a lower soft-tissue 
attenuation will be assumed in low weight and high height subjects. 

[0125] In another embodiment, a standard calibration curve is applied to x-ray 
10 images, whereby said calibration curve can be derived from reference x-rays obtained 
with use of calibration phantoms. For example, 100 patients can undergo dental x-rays 
with a calibration phantom and a standard calibration curve can be derived from these 
images. Similarly, 100 patients can undergo hip x-rays with a calibration phantom and a 
standard calibration curve can be derived from these images. Different calibration 
15 curves can be generated for different populations, for example, by generating different 
calibration curves for different ranges in body mass index, body height, sex, race etc. 
3.1.0.0. Anatomical Landmarks 

[0126] In one embodiment, identification of anatomic landmarks of the structure 
to be analyzed or identification of anatomical landmarks adjacent to the structure to be 

20 analyzed with subsequent positioning and computer analysis of the x-ray image relative 
to these anatomic landmarks or with subsequent positioning and computer analysis of 
anatomical region of interest (ROI) relative to these anatomic landmarks is performed. 
The present invention includes also positioning dental or other x-ray detectors, 
positioning the dental or other x-ray tube, and analyzing the resulting images using 

25 landmarks based on either 1) textural information, 2.) structural information, 3.) density 
information (e.g. density), or 4) 2 or 3 dimensional contour information 5) a 
combinations thereof of the tissue or structure to be measured and of tissues or 
structures adjacent to the measurement site. The invention also includes methods and 
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devices that are not necessarily based solely on anatomical landmarks, but in some 
applications can be combined with anatomical landmark embodiments. Preferably, 
many of the embodiments described herein are designed for automated use with a 
minimum of operator intervene and preferably remote or computer control of such 
5 devices. 

[0127] In one embodiment, an alignment device may be used to ensure 
perpendicular or near perpendicular alignment of the dental or other x-ray tube relative 
to the x-ray film, thereby decreasing geometric distortion resulting from tube angulation. 
For example, an x-ray film holder is positioned relative to an anatomical landmark, e.g. 

10 the posterior wall of the mandible in the incisor region. A side-view of an exemplary 
alignment system using a dental x-ray film holder is shown in Figure 1 5. The system 
includes bite block (100), stainless steel rod (101), film (103), optional calibration 
phantom (104), Rinn holder (102) typically having a ring or donut shape, and extension 
tubing (200). The extension tubing is designed to fit within the Rinn holder and may be 

15 temporarily or permanently attached. The system can achieve high reproducibility of 
the film position relative to an anatomical landmark such as the alveolar ridge or the 
posterior wall of the mandible. The extension tubing allows for alignment of the x-ray 
tube so that it is near perpendicular to the Rinn instrument and, ultimately, the dental 
film. 

20 [0128] Since manual alignment of the dental x-ray tube, namely the tube {e.g., 
metal) located in front of the dental x-ray tube for pointing and alignment purposes, is 
often not very accurate with alignment errors of 3, 5 or even more degrees, a 
mechanical or electromagnetic device is preferably used in order to achieve 
perpendicular or near perpendicular alignment between the metal tube anterior to the x- 

25 ray tube and the Rinn holder. For example, the metal tube can be physically attached 
to the Rinn holder with use of one or more Velcro™ straps or it can be aligned using 
optical aids such as levels, cross-hairs, light sources (points or areas), etc. 
Alternatively, such physical attachment can be achieved with use of one or more 
magnets rigidly attached to the dental x-ray system metal tube and the Rinn holder. In 
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this embodiment, the magnets on the Rinn holder and the dental x-ray system metal 
tube will be aligned and brought into physical contact. In another embodiment, an 
extension tube is attached, for example with an adhesive, to the Rinn holder. The 
extension tubing can also be an integral part of the Rinn holder. The extension tubing 

5 can be designed so that its inner diameter is slightly greater than the outer diameter of 
the dental x-ray system metal tube. The dental x-ray system metal tube is then inserted 
into the extension tubing attached to the Rinn holder thereby greatly reducing alignment 
error of the x-ray tube relative to the x-ray film. Alternatively, the extension tubing can 
be designed so that its outer diameter is slightly smaller than the inner diameter of the 

10 dental x-ray system metal tube. The dental x-ray system metal tube is then advanced 
over the extension tubing attached to the Rinn holder thereby greatly reducing 
alignment error of the x-ray tube relative to the x-ray film. One of skill in the art will 
easily recognize in view of the teachings herein that many other attachment means can 
be used for properly aligning the dental x-ray tube with the dental x-ray film. 

15 Combinations of attachment mechanisms are also possible. 

[0129] The anatomical landmark that is selected is part of an anatomical region. 
An anatomical region refers to a site on bone, tooth or other definable biomass that can 
be identified by an anatomical feature(s) or location. An anatomical region can include 
the biomass underlying the surface. Usually, such a region will be definable according 

20 to standard medical reference methodology, such as that found in Williams et al., Gray's 
Anatomy, 1980. The anatomical region can be selected from the group consisting of an 
edge of the mandible, an edge of the maxilla, an edge of a tooth, valleys or grooves in 
any of these structures or combinations thereof. The dental x-ray image can be readily 
taken so as to include the anatomical site. Other anatomical regions include but are not 

25 limited to the hip, the spine, the forearm, the foot, and the knee. 

[0130] For example, the region of interest is placed between the dental apices 
and the inferior mandibular cortex. The apices can be found automatically in the 
following way: for each row of pixels, the gray value profile is examined. While a profile 
that intersects bone and dental roots in an alternating fashion has several distinct peaks 
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and valleys, a profile that only covers trabecular bone shows irregular changes in the 
gray values (Fig. 6). The dental apices are located in the transitional region between 

these two patterns. 

[01 31] The measurement techniques to assess trabecular bone structure or 
5 macro-anatomical and/or biomechanical parameters are preferably designed to work 
without user intervention. In order to fully automate the process of analyzing dental x- 
rays, it is necessary to develop a technique to locate the regions of interest (ROIs) that 
are used for the calculation of the structural parameters of the trabecular bone. If the 
profile for a particular row of pixels contains distinct peaks, their number, width and 
10 height can be determined. Next, the rows below these lines can be evaluated until the 
peaks have disappeared. This line determines the boundary, 5 mm below which the 
ROI can be placed in the center between the longitudinal axes of the roots, which can 
also be determined from the row profiles (Fig. 6). At a pixel size of 0.042mm x 
0.042mm, which corresponds to a resolution of 600dpi, the ROI has a size of 5.4mm x 
15 5.4mm (128x128 pixels). For other scanning resolutions, the pixel resolution of the ROI 
can be adjusted accordingly. 

[01 32] In the case of an edentulous patient, bone mineral density can be 
measured in all ROIs that are located on a line that is, for example, 8 mm inferior and 
parallel to the alveolar ridge. The ROIs can be moved from left to right on a pixel-by- 

20 pixel basis. Eventually, the ROI with the lowest BMD can be chosen for further 

evaluation of the structural bone parameters. This helps to avoid inclusion of regions on 
the x-ray where bone mineral density may be overestimated due to projection of the 
curved parts of the mandible near the canine teeth. Alternatively, the ROI with the 
median BMD can be used. Other statistical parameters can be employed for this 

25 purpose. 

[01 33] Thus, software or other computational unit can identify the selected 
anatomic landmark in an interrogated x-ray image and direct analysis of the image 
using various parameters and analytic functions. Further, such software or other 
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computational analytical unit can be used to identify areas of particular density at a 
certain distance from the selected landmark. Similarly, manual or computer analysis 
can be used to identify areas of lowest, highest, median or average density (or 
structural charactsrtetics) in relation to the selected landmark. 

5 [01 34J Further, the same landmark may be compared at different times (intra- 
landmark comparison) or one or more landmarks may be compared (inter-landmark 
comparison). For instance, an intra-landmark comparison can be used during a single 
interrogation protocol that entails multiple interrogations of the same region with 
reference to a particular anatomical landmark. Statistical analysis as described herein 

10 and known in the art can be performed. 

[0135] Thus, the invention provides for means of assessing bone structure, i.e. 
the two-dimensional or three-dimensional architectural organization of the trabecular 
bone including, but not limited to, measurement of trabecular spacing, trabecular 
thickness, trabecular length and trabecular connectivity. Other examples of 

15 measurements of bone structure are provided in Table 1. These measurements can be 
used alone or enhanced with use of calibration phantoms or external standards that can 
allow a correction or normalization of image intensity and that can in certain 
embodiments also allow a correction of geometric distortions for example resulting from 
cone beam geometry of an x-ray beam. 

20 [0136] The invention provides for means of assessing macro-anatomical and/or 
biomechanical parameters. These measurements can be used alone or enhanced with 
use of calibration phantoms or external standards that can allow a correction or 
normalization of image intensity and that can in certain embodiments also allow a 
correction of geometric distortions including magnification, for example resulting from 

25 cone beam geometry of an x-ray beam. 

[01 37] As described herein, one or more measurements of bone structure or 
macro-anatomical and/or biomechanical parameters can be used to select a therapy, for 
example the use of anabolic or antiresorptive agent in the case of bone loss or 
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deterioration. In certain embodiments, measurements of bone structure and/or one or 
more macro-anatomical and/or biomechanical parameters are conducted overtime to 
longitudinally monitor a subject's bone health longitudinally over time. Measurements 
can be performed at different time points T1, T2, Tn and changes in said bone 

5 structure and/or macro-anatomical and/or biomechanical parameters can be registered 
and used to track a patient's bone health. In either single or longitudinally 
measurements, a physician can be apprised of the measurements and can include a 
pre-determined cut-off value (e.g., when a bone structure or macro-anatomical and/or 
biomechanical parameter measured in a patient is more than one or two standard 

10 deviations different from a normal, healthy reference population) and use this 
information to select a therapy. 

[0138] The data obtained and analyzed as described herein can be used to 
monitor a patient's response to therapy. For example, information regarding bone 
structural and/or macro-anatomical and/or biomechanical information in a patient 

15 receiving an anabolic or antiresorptive drug and be evaluated at different time intervals 

T1, T2 Tn and changes in said bone structure and/or macro-anatomical and/or 

biomechanical parameters can be used in order to assess therapeutic efficacy. A 
physician can use this information to adjust the dose of a drug administered (e.g., for 
treatment of osteoporosis) or to change the drug regimen. 

20 [0139] Other techniques using x-ray information such as tomosynthesis can also 
be used for measuring bone structure and for selecting said therapy or monitoring said 
therapy. 

[0140] Bone structure can be measured using a number of different technical 
approaches. These include but are not limited to the Hough Transform, analysis of 
25 density and size distribution of trabeculae, multidimensional classification schemes, 
mean pixel intensity, variance of pixel intensity, Fourier spectral analysis, fractal 
dimension and morphological parameters. 
3.1.1.0. Hough Transform 
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[0141] The Hough transform (See, e.g., Hough "Machine analysis of bubble 
chamber pictures" in International Conference on High Energy Accelerators and 
Instrumentation. 1959. CERN) can be used to detect geometric objects in binary 
images. As an entirely new approach to assessing bone structure and/or macro- 
5 anatomical, the invention includes the use of such methods to analyze direction and 
length of structures in bone images. For this purpose, the region of interest (ROI) can 
be blurred with a Gaussian filter. The pixel values of the filtered ROI can then be 
subtracted from those in the original ROI, and the value 128 can be added at each pixel 
location. This results in an image with a mean gray value of 128, which is also used as 
10 a threshold to yield a binary image in which the trabeculae are represented by the white 
pixels. 

[0142] After a skeletonization step, a Hough transform with the line 
parameterization p = xcosd + ysm6 can be applied to the binary image in order to find 
straight line segments. Here p is the perpendicular distance of the line from the origin 

15 and 6 is the angle between the x-axis and the normal. Each point (x,y) in the original 
image is transformed into a sinusoidal curve p = xcos0 + ysin9 in the (p,0) plane of the 
transformed image (see Fig. 7)). Ideally, the curves from collinear points in the original 
image intersect in a single point in the transformed image. However, the (p,&) plane 
can be divided into bins, where each bin counts the number of transformed curves that 

20 pass through it. This number corresponds to the number of collinear points on a line 
segment in the original image, and thus the length of this segment. Furthermore, the 
transformed image provides information on the predominant angles of the line segments 
in the original image (see Fig, 8). 

[0143] The average length and the variance of the line segments, which are 
25 calculated for all bins with a count above a certain threshold, can be used as structural 
parameters for the shape of the bone trabeculae. Average length as well as the 
variability of the length to decrease in patients with osteoporosis. The threshold has the 
effect that only segments of a certain minimal length are included in the calculation. 
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Choosing the threshold so that it provides the best discrimination between healthy and 
diseased individuals can be readily determined by one of skill in the art in view of the 
teachings herein. 



in which each bin is interpreted as an element with a mass equivalent to its count, is a 
way to measure the predominant angles of the trabecular segments. The angle at cm is 
measured with respect to the alveolar rim to obtain a standardized value. More 
importantly, the variance of the segment angles (again measured after thresholding the 

10 bin counts) provides information on the anisotropy of the trabecular structure. 

Histomorphological studies of osteoporotic vertebrae have shown that the variability of 

trabecular orientations decreases with the disease. 

3.1.2.0. Analysis of Density and Size Distribution of Trabeculae 

[0145] Morphological operations such as variations of dilation and erosion and 

15 combinations thereof can also be used to detect the size of structures in gray scale or 
binary images. For example, a skeleton operator can be used to extract and quantify 
trabeculae of different sizes and directions, which results in a measure of the size 
distribution of trabecular structures. This skeleton operator is based on the work 
described in Kumasaka et al. (1997) Dentomaxillofac Rad 26: 161 -168 and works as 

20 follows: 

[0146] Let a two-dimensional structuring element e be a function over the window 
-m<, i, j < m (m>0) with E(i, j) e {o,l} . The dilation operator sets a pixel value f(x,y) in a 
gray scale image f to the maximum of those values within the window of size m, for 
which e(/j)=1: 



The "center of mass" of the transformed image h, given as, 



5 
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[f © E\x, y) = J%*Jf(x + i,y + J)\W, j) = l} 

[0147] The erosion operator is defined accordingly, using the minimum instead of 
the maximum: 

\f ® E\x, y) = _ min Jf(x + i,y + j)\E(i, j) = l} 
5 [0148] 'Opening' is the operation of maximum search after minimum search: 

f E =(f®E)®E 

[0149] Accordingly, the 'closing' operation is defined as the minimum search after 
maximum search: 

f E ={f@E)®E 

10 [0150] If a fixed structuring element E1 is given as E 1 (/j)=1 for -l<i,;<l, the 
skeleton operation is then defined as 

Sw— (f) = (/ ® Ei ) - [(/ ® £2)3 ] (V 

[0151] E 2 is another structuring element that is of circular shape and can be 
15 varied in size, and therefore renders the skeleton operator sensitive to the size of the 
structures in the image. The erosion of f with E 2 erases the structures that are smaller 
than E 2 and extracts those trabeculae that are at least equal in size. Those structures 
that are exactly equal in size is reduced to a width of one pixel. The opening step with 
E-i causes all structures that are one pixel wide to disappear (second term in (1)). After 
20 subtraction of this term from the first one, only those trabecular structures that exactly 
match the size of E 2 remain. Finally, the image is thresholded with a level of 1 . The 
effect of this operator is illustrated in Fig. 9. 

[0152] Fig. 10 demonstrates the use of the skeleton operator with the same 
structural element diameters as in Fig. 9 on a gray scale region of interest from a dental 
25 x-ray containing trabecular bone. The number of bright pixels in the binary images 
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resulting from each skeleton operation corresponds to the portion of trabeculae of the 
particular size in the original image. If the percentage of the bright pixels with respect to 
the total number of pixels in each skeletonized image is plotted against the diameter of 

E 2 , the "center of mass" of the curve, i.e. the predominant structure size, can be used as 
5 an index to discriminate between osteoporotic and healthy bone. 

[0153] Furthermore, the skeleton operator is preferably optimized and extended 
to detect structures that are oriented only in a specific direction. This can be achieved 
by adding erosion operations to the skeleton operator with structural elements in which, 
for example, only the diagonal pixels are set to 1 . 
10 [0154] This can be used to calculate an anisotropy index, similar to the one 
derived from the Hough transform. Both anisotropy indices are tested with respect to 
their potential to distinguish healthy from osteoporotic bone. 

[0155] In a similar manner the sizes of the marrow spaces can be examined. The 
skeleton operator is then defined as 



subtracted gray level structures on x-ray images to characterize the homogeneity of 
trabecular structures. This process takes into account the gray level contrast between 
structures to define marrow spaces. The watershed segmentation, when applied to 

20 background subtracted bone x-ray images, defines regions with lower gray levels (or 
basins) surrounded by higher gray level structures (or ridges), as marrow space, in 
accordance to the spatial extend and gray levels of ridges. Therefore, the size and 
orientation of marrow space segments defined by this procedure can be related to the 
spacing, relative density and orientation of adjacent trabecular structures. The 

25 segments of marrow space generated using the watershed segmentation can be 

measured for their area, eccentricity, orientation, and the average gray level on the x- 
ray image within the segment. The statistics (for example mean, standard deviation, 



15 




[0156] 



In addition, the watershed segmentation can be applied to background 
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minimum, maximum, and mode) for each of the segment characteristics can be 
measured. These statistics can be selected to reflect the homogeneity of marrow space 
and trabecular structures, and can be used to detect presence of abnormal distribution 
of marrow space and trabecular structures. 
3.1.3.0. Multidimensional Classification Schemes 

[0157] In certain embodiments, it is preferred to use multiple indices to measure 
bone structure and/or macro-anatomical parameters. Thus, novel approaches that 
integrate one or more suitable indices can be employed. The indices can be optimized 
and incorporated into a multi-dimensional classification scheme, for example using a 
nearest neighbor classification. Cover et al. (1967) IEEE Trans Inform Theory 13(1):21- 
7. (See, Example 3). 

[0158] Table 1 provides examples of different analyses and anatomical / 
physiological correlates of the parameters that can be measured. 



Table 1 



Analysis 


Anatomical / Physiological Correlates 


Hough transform 


length and direction of trabecuiae; anisotropy 


Morphological 
operators 


thickness and direction of trabecuiae; anisotropy; thickness and 
length of marrow spaces 


Mean pixel intensity 


bone mineral density 


Variance of pixel 
intensity 


complexity of trabecular structure 


Fourier spectral 
analysis 


complexity of trabecular structure 


Fractal dimension 


complexity of trabecular structure 


Morphological 
parameters 


length, size of trabecuiae; complexity of trabecular structure; length, 
size of marrow spaces; complexity of marrow space 



15 3.1.3.1 MEAN PIXEL INTENSITY 

[01 59] Mean pixel intensity is a general parameter for the bone mineral density. 
The degree to which x-rays passing through bone tissue are absorbed depends on the 
bone's mineral content. Bone with a higher mineral density absorbs a larger portion of 
x-rays, and therefore appears brighter on an x-ray image. 
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[0160] The mean pixel intensity f(x,y) in the ROI is calibrated against an 
aluminum calibration wedge that is included in the image. The log of the average pixel 
intensity for each thickness level of the calibration wedge is plotted against the 
thickness, which allows f(x,y) to be converted into a standardized aluminum thickness 
5 equivalent, which is used as the value for this parameter. The automatic recognition of 
the different thickness levels of the calibration wedge are made possible by different 
geometric patterns scribed into the wedge which are shown in the x-ray image and can 
be localized automatically. 

3.1.3.2. Variance of pixel intensity 

10 [0161] The variance of the pixel gray values in the ROI, varf(x,y), describes the 
variability of the pixel intensities and can therefore be a measure of the degree of 
trabeculation. A loss of trabecular bone is predicted to be reflected by a decreased var 
f(x,y). Southard & Southard (1 992) Oral Surg Oral Med Oral Pathol 74: 1 1 1 -1 1 7. 

3.1.3.3. Fourier spectral analysis 

15 [0162] The spatial frequency spectrum of a texture provides information about its 
coarseness. Fine textural structures and edges in an image correspond to high 
frequencies in the frequency domain, while coarse textures are represented by lower 
frequencies. Applied to x-ray images of trabecular bone, this means that a region with 
coarse or little trabeculation should exhibit a Fourier spectral energy concentration at 

20 low spatial frequencies, whereas a region of fine trabecular structure should show a 
spectral energy concentration at high frequencies. 

[0163] Typically, the 2-dimensional Fourier coefficients for the selected ROI. 
These 2-dimensional coefficients are used to determine a 1 -dimensional power 
spectrum F(u) by averaging all coefficients over circles with radii that correspond to the 
25 discrete spatial frequencies u. The mean transform coefficient absolute value \F(u)\ and 
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±\F(u)\.u 

the mean spatial first moment M, = *~ 2 N l — of the absolute coefficients are 

determined after exclusion of the first ("DC") coefficient. Mi provides a measure for 
which frequencies contribute most to the energy of the spectrum, similar to the "center 
of mass" of a geometric object. 

5 

3.1.3.4. Fractal dimension 

[0164] A different approach to analyze the texture in an image is by fractal 
analysis. Fractals are objects that exhibit certain statistical self-similar or self-affine 
10 properties, so that a portion of the object, scaled to its original size, has for example the 
same surface area (3-d) or the same perimeter (2-d) as the original object. In the 
context of fractal analysis, the gray values in a particular texture can be interpreted as 
an altitude, and the resulting 3-dimensional surface is analyzed (Fig. 11). 

[0165] Fractal dimension (fd) is the rate at which the perimeter or surface area of 
15 an object increases as the measurement scale is reduced. Russ "The Image 

Processing Handbook," Third edition ed. 1999, Boca Raton: CRC press. It is a measure 
for the complexity of a boundary or surface and corresponds to the intuitive notion of an 
object's roughness. Without being bound by one theory, it is postulated that 
osteoporotic trabecular bone, in which trabeculae become thinner and lose their 
20 continuity, and therefore complexity is increased, should have a higher fractal 
dimension than healthy bone. 

[0166] The results from the several ways in which FD can be measured are not 
comparable. Thus, various methods can be tested to determine which one (or 
combination) provides the best discrimination between normal and osteoporotic 
25 subjects. 
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[01 67] The first method is applied in the frequency domain after calculation of the 
ROI's 2-D power spectrum using a fast Fourier transform (FFT). From the 2-D Fourier 
coefficients the 1-D power spectrum is produced as described above for the Fourier 
analysis. When this 1-D power spectrum is plotted as the logarithm of the power versus 
5 the logarithm of the frequency, it must have a negative slope of magnitude b with 1 < b 
< 3 according to fractal theory. The FD value is then calculated as fd 1 =3.5-6/2 . 

[0168] Another approach, the Minkowski method, measures the difference 
(summed over the ROI) between an upper and lower envelope fitted to the surface as a 
function of the size of the neighborhood used. Peleg et al. (1 984) Anal Mach Intel! 
10 6(4):518-523. If 8 (8=1,2,3,...) is the distance between the envelopes and the surface, 
then the upper envelope u 5 and the lower envelope i 8 are given by 

u 0 (i,j) = l 0 (i,j) = f(ij) 
« s+ i ft j) = max ja, ft j) + 1, ^ max ^ {u s (m, w)}J 
l M ft J) = minji, ft j) - 1, ^ {/, (m, «)}} 

where f(i,j) is the gray value of pixel (i,j) in the ROI. The log of the area A(8), plotted 
against log(8), yields a line with a negative slope of magnitude b'. The fractal dimension 
15 is then given by FD 2 -2-b' . The area is calculated as A{5) = Vs * s ~ ] with 
v,= D(« a ftj)-^ft/))- 

(/J)eB07 

3.1.3.5. Morphological Parameters 

[0169] While the previous features and parameters provide rather general 
information on trabecular bone structure, the following examples describe more detailed 
20 aspects. 

[0170] The gray scale region of interest is first binarized. As described in White 
et al. (1999) Oral Surg Oral Med Oral Patholo Oral Radiol Endod 88:628-635, this can 
be achieved in the following way: The ROI is blurred by means of a Gaussian filter. The 
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blurred ROI is then subtracted from the original ROI, and the value 128 is added at each 
pixel location. This results in an image with a mean gray value of 128, which is also 
used as a threshold, resulting in an image, in which trabeculae are white and marrow 

space is black. 

5 [0171] From this binary image, the total number of white pixels represents the 
trabecular area, which is calculated as a percentage of the total ROI area. The number 
of pixels on the outer trabecular border measures the peripheral length of the 
trabeculae. The same parameters can be measured for the marrow space by counting 
the black pixels. 

10 [0172] After skeletonization of the binary image, the total length of the trabeculae 
is determined by the total number white pixels. Furthermore, counts of the terminal 
points and of the branch points are expressed as a proportion of trabecular length. An 
estimate of the average length of the trabeculae is calculated as the ratio of total 
trabecular length and the sum of terminal points and branch points. 

15 3.1.3.5. Markov Random Fields 

[0173] In certain embodiments, Markov random fields can be used as models for 
osteoporosis detection from radiographic images and for fracture risk prediction. As 
noted herein, osteoporosis is typically manifested in radiographic images by structural 
changes that can be used for computer-aided detection and characterization. Thus, the 

20 detection and/or characterization of osteoporosis from radiographic images relies on the 
measurement and analysis of a feature or set of features relating to the density of the 
bone or trabecular structures present in an image. 

[0174] Markov random fields can be used to analyze and detect structure density 

changes by either modeling particular joint feature distributions ({F1, F2 Fn}) as 

25 they are estimated at each image element or image neighborhood (Figure 1 9), or by 
modeling the actual radiographic manifestation of particular structural definitions 
(e.g. trabeculae) (Figure 20). 
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[0175] In the first case of estimation at each image element or image 
neighborhood, the Markov random field framework is used for a context-based feature 
analysis/discrimination approach which takes into account local relationships between 
the features and effectively compensating for space-varying processes (e.g. variable 
5 soft tissue or missing or incomplete data due to boundaries) that can affect the relative 
values of the features taken into account. (Buntine (1994) "Operations for learning with 
graphical models," J. Artificial Intelligence Res. December:! 59-225). 
[01 7@] This approach can also be used for predicting most likely fracture paths 
based on the analysis of trabecular structure nodes and their related feature sets by 
10 defining the most likely chains of joint feature sets. The analysis framework can be a 
^ _ P(ROI | Normal) 
Likehood Ratio approach: p ( ROI | Abnormal )_ wnere P( RO | | . ) is given by 

the corresponding Markov random field model. 

[0177] Another analysis approach is through the implementation and training of 
Bayesian networks, for example as described in Heckerman D (1996) "A tutorial on 
15 learning with Bayesian networks," Microsoft Research Technical Report, MSR-TR-95- 
06, based on available test case data. 

[0178] Markov random fields can also be used to model the manifestations of the 
structures in an image in probabilistic terms. (Geman et al. (1984) "Stochastic 
relaxation, Gibbs distributions, and the Bayesian restoration of images," IEEE 

20 Transactions on Pattern Analysis and Machine Intelligence 6:721 -741 ; Besag (1 986) 
"On the statistical analysis of dirty pictures," Journal of the Royal Statistical Society, 
48(3):259-302). As depicted in Figure 20, each of the image components (noise and 
characteristic structure) has associated probabilistic models, P(N) and P(l | T) 
respectively, that describe the spatial distribution of the gray-level intensity. For 

25 example, a common assumption for the noise component in digital/digitized radiographs 
is to consider Normal or Poisson distributed pixels. The nomenclature for the 
distribution of the characteristic texture P(l | T) is such as to reflect that the 
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corresponding probability distribution of the region I is conditional ( expressed by the 
symbol | )on the characteristic structure present T. The analysis tools for such a 
probabilistic framework are provided by the laws of probability and specifically Bayes' 
Rule shown in Figure 21 . Bayes' rule can be described as the rule according to which 
5 our knowledge about the presence of a given characteristic structure in an ROI is 

updated (a-posteriori information represented by the probability distribution P(T | ROI) ), 
based on experience of how often (or likely) each characteristic structure is present (a 
priori information represented by the probability distribution P(T) ) and knowledge of 
how the sources of noise and variability change the manifestation of the corresponding 

10 characteristic structure (knowledge of the likelihood, thus also called the likelihood 
function P( ROI | T ), of the ROI image given the possible characteristic structures and 
overlapping degrading components). Figure 21 illustrates that simply selecting the 
structure with the maximum a-posteriori information can be used as a decision criterion. 
[0179] To define the likelihood function P(ROI | T), Markov random field modeling 

15 may be employed. Markov random fields are specific multidimensional random 
processes that satisfy what is known as the Markov property. The Markov property 
simply states that in a random series of events, each event can be predicted and 
depends only on a limited set of events. This property is convenient and intuitive for the 
modeling and analysis of structures in images. It basically states that if the distribution 

20 of pixels in an ROI can be modeled as having the Markov property, then in order to 
determine if a pixel belongs to a given structure, only a limited number of neighboring 
pixels are necessary. 

[0180] Random fields having the Markov property confer the additional benefit of 
having an associated Gibbs probability distribution given by the following equation: 

25 

-U( Sl) s 2I ...) 
P( ROI = s 1)S2I ... ,sj=-| 
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[0181] Where the function: 



U(s 1 ,s 2I ...) = 2V(s i ) + 2V(s.,s.) 

+ 2V(s iI s j .,s k ) + ... 

5 depends on functions V (called potentials) of local neighboring elements called cliques: 
{s i },{s i ,s.},{s |I s J) s h } > ... e C 

[0182] The significance of cliques is that they are the fundamental elements that 
can be used to reflect specific spatial distribution properties of a structure of interest, 
10 such as for example vertical, horizontal and diagonal geometries. Furthermore, the 
Markov property is manifested very conveniently as each image pixel can be expressed 
in terms of the cliques in a local neighborhood: 



e V - 

P(s,|s k ,s k£ N,) = 



s c 

[01 83] The model parameterization for the families of images characteristic of a 
particular structural density grade and definition of a priori information can be done 
either by estimation from available patient data thus defining empirical priors or by 
implementing physical and stochastic models that are based on the image generation 
process. 

3.1.4.0. Overlapping windows processing 

[0184] Furthermore, two or more overlapping ROIs can also be defined and used 
to analyze any given image. In other words, bone density, microarchitecture , macro- 
anatomic and/or biomechanical (e.g. derived using finite element modeling) analyses 
can be applied within a region of predefined size and shape and position. This region of 
interest may also be referred to as a "window." Processing can be applied repeatedly 
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within the window at different positions of the image. For example, a field of sampling 
points may be generated and the analysis performed at these points (Fig. 16). The 
results of the analyses for each parameter can be stored in a matrix space, e.g., where 
their position corresponds, to the position of the sampling point where the analysis 
5 occurred, thereby forming a map of the spatial distribution of the parameter (a 

parameter map). The sampling field can have regular intervals or irregular intervals with 
varying density across the image. 

[0185] The amount of overlap between the windows can be determined, for 
example, using the interval or density of the sampling points (and resolution of the 

10 parameter maps). Thus, the density of sampling points is set higher in regions where 
higher resolution is desired and set lower where moderate resolution is sufficient, in 
order to improve processing efficiency. The size and shape of the window would 
determine the local specificity of the parameter. Window size is preferably set such that 
it encloses most of the structure being measured. Oversized windows are generally 

15 avoided to help ensure that local specificity is not lost. 

[0186] The shape of the window can be varied to have the same orientation 
and/or geometry of the local structure being measured to minimize the amount of 
structure clipping and to maximize local specificity. Thus, both 2D and/or 3D windows 
may be used, depending on the nature of the image and data to be acquired. 

20 [0187] In another embodiment, bone density, microarchitecture, macro-anatomic 
and/or biomechanical (e.g. derived using finite element modeling) analyses can be 
applied within a region of predefined size and shape and position. The region is 
generally selected to include most or all of the anatomic region under investigation and, 
preferably, the parameters can be assessed on a pixel-by-pixei basis (e.g., in the case 

25 of 2D or 3D images) or a voxel-by-voxel basis in the case of cross-sectional or 

volumetric images {e.g., 3D images obtained using MR and/or CT). Alternatively, the 
analysis can be applied to clusters of pixels or voxels wherein the size of the clusters is 
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typically selected to represent a compromise between spatial resolution and processing 
speed. Each type of analysis may yield a parameter map. 

[0188] Parameter maps can be based on measurement of one or more 
parameters in the image or window; however, parameter maps can also be derived 
5 using statistical methods. In one embodiment, such statistical comparisons can include 
comparison of data to a reference population, e.g. using a z-score or a T-score. Thus, 
parameter maps can include a display of z-scores orT-scores. 
3.1.4.1. Analysis and selection of parameter maps 

[0189] The parameter maps can represent individual parameters or combinations 

10 of parameters such as density, microarchitecture macro-anatomical parameters or 
biomechanical parameters, for example derived using finite element modeling, are 
useful in identifying regions or patches that have similar characteristics. For instance, 
depending on their position, shape, size, orientation, and extent particular regions or 
patches that exhibit similar characteristics (e.g., values at high or low ranges of the data 

15 set) typically represent regions of bone with different properties, for example areas of 
stronger or weaker areas. Therefore, parameter maps can be used to generate virtual 
fracture lines that aid in predicting areas of the bone that may be subject to an 
increased risk of fracture. One or more parameter maps can be selected by statistical 
analysis of results from in vitro mechanical loading tests or by other means (e.g. from 

20 cross-sectional or longitudinal studies in osteoporosis subjects, in particular those 

developing fractures). Selection can be based, for example, on patch location, shape, 
size, orientation and extent that best correlates with location of actual fracture lines 
and/or for having parameter values that are best correlated with fracture risk, the 
incidence of osteoporotic fractures or fracture loads. 

25 3.1.4.2. Fracture path prediction 

[01 90] When there are multiple parameter maps that correlate well with fracture 
line, a multivariate regression model can be fitted to generate a composite parameter 
map derived from 2D or 3D data sets, e.g. x-rays, digital tomosynthesis, CT and MRI, 
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using the techniques described herein and/or statistical methods known to those of skill 
in the art. A parameter map can be used to predict the overall bone strength or fracture 
risk or fracture load by analyzing the predicted fracture paths. A predicted fracture path 
is defined here as the hypothetical path where fracture would most likely to occur, if 

5 sufficient forces are applied in one or more particular directions. 

[0191] In certain embodiments, a watershed segmentation can be applied to the 
selected or composite parameter map. Watershed segmentation can be applied to 2D 
images as well as to 3D (cross-sectional or volumetric data obtained, for example, from 
CT or MR). The boundaries of watershed segmentation generally form along the 

10 ridges on the parameter map, i.e., along the peak values. For a parameter that is 

positively correlated to bone strength or fracture load, i.e., higher values correspond to 
stronger bones, the inverse value of the parameter is used to generate the watershed 
boundaries so that the boundaries would form along valleys (local minimum) of 
parameter maps. The nodes of watershed boundaries can be identified and segmented 

15 to separate the watershed boundaries into segments (Figure 17). Each of these 

segments can be assigned a strength value or fracture load value which is a composite 
value of one or more parameter maps underlying the segment. The length, orientation, 
and position of segments can be used as normalizing factors for the strength values. 
[01 92] The nodes and segments of the watershed boundaries may be labeled, 

20 traced, measured, and recorded in a form of data structure, for example, a graphical 
structure. The strength values and interconnect relationships are also stored for each 
segment. To identify the most likely fracture paths, a search strategy, for example, the 
depth-first search (Russell S., Norvig, P., Artificial Intelligence: A modern approach. 
1995, NJ: Prentice Hall, pp.77), is propagated through the data structure to determine 

25 the paths of least resistance from one surface of the bone to another opposite surface 
restricted by a predefined solid angle. Alternatively, an artificial neural network can be 
trained to predict fracture paths given the parameter maps as inputs. 
3.1.4.3. Fracture risk prediction 
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[0193] Having predicted one or more fracture paths, additional processing may 
be performed, typically with a new processing grid that has high concentration of nodes 
along the predicted fracture paths with a different window size and/or shape. Macro- 
anatomical parameters such as cortical thickness can be evaluated (in two or three 
5 dimensional images) with higher resolution at the exits of fracture paths. Parameters 
that are the best predictors of fracture risk can be evaluated along the predicted fracture 
paths. These parameters, including density, microarchitecture, macro-anatomical 
measurements and biomechanical parameters, are selected by statistical analysis of 
results from in-vitro mechanical loading test or by other means, e.g. using cross- 
ly sectional or longitudinal studies in osteoporosis subjects, in particular those developing 
fractures, for being highly correlated to the magnitude of one or more mechanical 
properties of bone, for example in one or more particular loading force directions, or for 
being highly correlated with fracture risk, incidence of new fractures or fracture loads. 
The mechanical properties include but are not limited to yielding load, stiffness, and 
15 Young's modulus. 

[01 94] The values of parameters along the predicted fracture paths may be 
compared against the statistical distribution of the population. The z-score and T-score 
of each parameter relates to the risk of fracture occurring in a particular predicted 
fracture path. Thus, a fracture risk score can be assigned to that fracture path. The 
20 predicted fracture paths can also be associated with the clinical definition of common 
fracture types. The overall fracture risk can then be evaluated by weighing fracture risk 
score of each predicted fracture path with the probability of a particular type of fracture 
occurring. Figure 18 depicts an exemplary summary of this process. 
3.1.5.0. Biomechanical assessment 
25 [0195] The features and values extracted from the processing of density, micro- 
architecture, macro-anatomical parameters can be used as the inputs for biomechanical 
modeling, for instance modeling using finite element analysis. Finite element modeling 
(FEM) can be used as a surrogate for the physical mechanical properties of bone or 
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composite of bone and implants. Briefly, FEM involves the division of a structure or 
object into discrete shaped elements, where the mechanical behavior of each element 
can be described by precise mathematical equations. Structural finite element analysis 
(FEA), e particular subset of FEM, is the calculation of the mechanical behavior (stress 
5 and strain) at any point within the structure under specific loading conditions. The 
foundation of every finite element model is the two-dimensional or three-dimensional 
data of the object or structure 

[01 96] Examples of microarchitecture and micro-anatomical features that can be 
used as input mesh for finite element analysis include but are not limited to the actual 

10 and derivation of image or data structures of trabecular structures, image or data 

structures of cortical bone, image, data structures of trabecular skeleton or parameter 
maps derived from overlapping window processing. As described herein, the input 
features can be obtained from 2D and/or 3D images. The application of simulated force 
can be in one or more directions, and is typically associated with the actual force 

15 components that would occur in a fracture incident. The finite element analysis 

provides an estimate of load and direction of fracture for each fracture incident scenario. 
Fracture risk is estimated by weighing the fracture loads with the probability of each 
fracture scenario occurring. Further, the fracture paths estimated by finite element 
analysis can be used as inputs to the analysis of density, micro-architecture, macro- 

20 anatomical features. For example, density, micro-architecture, macro-anatomical 
features can be measured in areas of fracture paths predicted by finite element 
modeling. Conversely, finite element analysis can be combined with additional image 
and clinical data to determine fracture risk by predicting if the bone would fracture, given 
the force components that would occur in a fracture incident. 

25 [0197] Bone fracture risk can be evaluated using one or a composite of more 
than one dependent or independent results of analysis or statistical methods. An 
example of this combination is the weighted average score of density, micro- 
architecture, macro-anatomical, finite element analysis and clinical risks factors such as 
weight, height, history of fracture, family history of fracture, and the like. 
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[0198] Finite element modeling can be applied to all of the bony structures 
included in an image. Preferably, however, finite element modeling is typically applied in 
selected subregions. In certain embodiments, finite element modeling is applied in 
areas coinciding with or bordering with the predicted fracture path, for example based 

5 on micro-structural or macro-anatomical measurements. By combining biomechanical 
assessment of bone properties with density, micro-architectural and macro-anatomical 
assessment, the prediction of fracture risk and/or the correlation with fracture load can 
be improved. Finally, regional assessment of biomechanical properties can also 
improve the accuracy of the fracture path prediction. 

10 [0199] Biomechanical assessment can also include more traditional approaches 
estimating levers and forces at the macro-anatomical level, e.g. measurement of 
moments, shear and compressive forces based on macro-geometric parameters of the 
bone and anticipated loads or stresses. These more traditional approaches can be 
combined with finite element modeling, measurements of density, bone structure, and 

15 macroanatomical parameters, e.g. cortical thickness, thereby improving assessment of 
bone strength and fracture risk and improving the correlation with fracture loads and, 
ultimately, incident new fractures. 

[0200] As will be appreciated by those of skill in the art, the macroanatomical 
parameters that are measured can change depending on the region of interest to be 

20 measured. For example, when studying a portion of the spine, the user can combine 
bone structure measurements with macroanatomical measurements and/or FEA and/or 
other biomechanical measurements and/or bone mineral density. The actual 
macroanatomical measurements that are used in the spine can be, for example, the 
inner pedicle distance, the outer pedicle distance, the vertebral height (either anterior, 

25 central, posterior, left, right, or a combination thereof), the vertebral anterior-posterior 
diameter (taken either in the superior, middle, inferior, or another location), the vertebral 
right to left diameter (taken in either the superior, middle, inferior or another location), 
the vertebral diameter (taken in an oblique plane), the vertebral diagonal (using, e.g., 
internal cortex or external cortex), the thickness of the superior endplate (taken, e.g., 
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anteriorly, centrally, posteriorly, from the left, from the right, or a combination thereof), or 
using the thickness of the inferior endplate (again taken, e.g., anteriorly, posteriorly, 
from the left, from the right, or a combination thereof). 

[0201] Similarly, when studying the knee and tibia, the user can combine bone 
5 structure measurements with macroanatomicai measurements and/or FEA and/or other 
biomechanical measurements and/or bone mineral density. However, as will be 
appreciated by those of skill in the art, the bone structures used for measurements 
when studying the knee and tibia region change due to changes in anatomy. Thus, in 
studying the knee and tibia region, suitable measurements are taken from, for example, 

10 the anterior-posterior diameter of the bone using the inner or outer cortex, or a 

combination thereof, the medial-lateral diameter of the bone using the inner or outer 
cortex, or a combination thereof, the cortical thickness in various locations, the standard 
deviation of cortical thickness, the subchondral bone thickness in various locations, 
and/or a combination thereof. 

15 [0202] Cases may arise where the macroanatomicai measurements are used to 
normalize bone structure or bone density measurements. For example, in the tibia, 
bone structure and/or bone density measurements could be altered if the patient has a 
tibia that is thick in the anterior-posterior dimension (e.g., thicker than average). The 
macroanatomicai measurements are then used to normalize the tibial measurement by, 

20 for example, forming a ratio between the thick tibial measurement in the anterior- 
posterior direction and another measurement. 
3.2.0.0. Soft Tissue 

[0203] Variations in soft tissue thickness can be significant in analyzing and 
evaluating bone density and bone structure, macro-anatomical parameters and 
25 biomechanical parameters, e.g. those derived using finite element modeling, in x-rays. 
Accordingly, the invention also includes methods and devices for correcting for soft 
tissue in assessment of bone structure or dense tissue, particularly for diagnosing 
and/or predicting osteoporosis or other bone conditions. 



60 



WO 2004/086972 



PCT7US2004/00<M<>5 



[0204] In certain embodiments, the x-ray image is a dental x-ray image and such 
correction methods involve (a) interrogating at least a portion of a subject's mandible 
and/or maxilla with an x-ray detector; (b) producing an x-ray image of the interrogated 
mandible and/or maxilla; (c) obtaining data from the x-ray image regarding bone density 

5 or bone structure; (d) interrogating the surrounding soft tissue to determine soft tissue 
thickness; and (e) correcting the data obtained from the x-ray image by correcting for 
soft tissue thickness. Such study groups include: non-osteoporotic premenopausal, 
non-osteoporotic postmenopausal, osteoporotic postmenopausal patients. It will be 
apparent, although exemplified with respect to dental x-rays, that many of the methods 

10 described herein can be applied to other x-ray images, e.g. hip or spine x-ray images. 
[0205] Soft tissue thickness measured in a subject can also be compared to 
reference soft tissue thickness obtained from a control population (e.g. age-, sex-, race-, 
or weight-matched normal subjects). Reference soft tissue thickness can be generated 
by measuring soft tissue thickness in healthy subjects with normal vascular, cardiac, 

15 hepatic, or renal function and no other underlying medical condition. Reference soft 
tissue thickness can be expressed as but are not limited to, mean and standard 
deviation or standard error. Reference soft tissue thickness can be obtained 
independently for patients 15-20, 20-30, 30-40, 40-50, 50-60, 60-70, 70-80, and 80 and 
more years of age and are preferably obtained separately for men and women and for 

20 race (e.g. Asian, African, Caucasian, and Hispanic subjects). Additionally, reference 
soft tissue thickness can be obtained for different subject weights within each age, sex, 
and racial subgroup. 

[0206] Individual patients can be compared to reference soft tissue thickness. If 
patient's soft tissue thickness is elevated, a correction factor can be applied. The 
25 amount/magnitude of correction factor is influenced by the magnitude of increase in soft 
tissue thickness that can be influenced by the magnitude of fat, fibrous, and muscle 
tissue contribution. Clinical study groups can be evaluated to generate databases for 
further study or to generate more refined correction factors. Such study groups include: 
non-edematous non-osteoporotic premenopausal, non-edematous non-osteoporotic 
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postmenopausal, non-edematous osteoporotic postmenopausal; edematous non- 
osteoporotic premenopausal, edematous non-osteoporotic postmenopausal, and 
edematous osteoporotic postmenopausal patients. In each study group the following 
procedures can be performed for comparison: dual x-ray absorptiometry ("DXA") of the 

5 spine, hip, or calcaneus, along with SOS and BUA measurements or quantitative 
computed tomography ("QCT"). Thus, correction for soft tissue thickness can also 
improve the accuracy and discriminatory power in the analysis of x-rays and other x- 
rays. Such methods can also be used to identify population with an increased or 
decreased risk of bone conditions such as osteoporosis. 

10 4.0. Applications 

[0207] The measurements of bone mineral density or trabecular architecture 
and/or macro-anatomical and/or biomechanical parameters, for example in the 
mandible or maxilla or in the hip or in the spine, can be used to derive an assessment of 
bone health in any subject. Additionally, the analysis and manipulation of data from x- 

15 rays allows for the assessment of bone health that in turn can be used to prescribe a 
suitable treatment regime. Efficacy of a treatment regime can also be assessed using 
the methods and devices described herein (for example, using measurements of bone 
mineral density or trabecular architecture and/or macro-anatomical and/or 
biomechanical parameters in the mandible or the maxilla or the hip or the spine taken at 

20 two separate time points T1 and T2 to detect any difference in bone mineral density or 
trabecular architecture). 

[0208] In addition, the methods described herein permit, for example, fully 
automated assessment of the structural organization and architectural arrangement of 
trabecular bone and/or macro-anatomical and/or biomechanical parameters on standard 
25 hip radiographs as well as improved tools for monitoring progression of osteoporosis 
and therapeutic response. In certain embodiments, the methods involve binarizing and 
skeletonizing trabecular bone using morphological operators with detection of branch 
points and endpoints of the skeleton network and classification into free-end segments 
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and node-to-node segments. In other embodiments, the methods involve measuring 
trabecular density, trabecular perimeter, trabecular bone pattern factor, segment count, 
segment length, angle of segment orientation and ratio of node-to-node segments to 
free-end segments based on the binarized and/or skeletonized images. In still further 

5 embodiments, the methods involve (a) measuring trabecular thickness using a 
Euclidean distance transform (see, also Example 3); (b) assessing trabecular 
orientation using a 2D Fast Fourier Transform; and/or (c) creating a bone structure 
index for diagnosing osteoporosis or for predicting fracture risk combining at least two or 
more of these structural parameters. 

10 [0209] In certain embodiments, the radiograph is of a subject's hip. Furthermore, 
to help control the influence of radiographic positioning on the accuracy of bone 
structure and/or macro-anatomical and/or biomechanical measurements, the methods 
may include one or more of the following: evaluating the angular dependence of bone 
structure measurements in the hip, for example by comparing antero-posterior 

15 radiographs of the hip joint in healthy to osteoporotic patients (subjects) with the femur 
radiographs in neutral position and in various degrees of internal and external rotation 
or by obtaining radiographs of the hip with different degrees of tube angulation. Bone 
structure and/or macro-anatomical and/or biomechanical measurements can be 
compared between the different positions to determine which bone structure parameters 

20 show the least dependence on radiographic positioning and/or using a foot holder to fix 
the patients' foot in neutral position in case pair wise coefficients of variation between 
the results for the 0° neutral position and a 15° internal or external rotation position 
exceed 10% for the majority of the structural parameters measured. 
[0210] In other embodiments, methods of monitoring bone structure and/or 

25 macro-anatomical and/or biomechanical parameters over time {e.g., longitudinally) are 
also provided, for example to assess progression of osteoporosis and/or response to 
therapy. In certain embodiments, the methods involve automated placement of regions 
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of interest (ROI) in the hip joint, for example by creating and using a general model of 
the proximal femur that includes six defined regions of interest (ROI's). 
[021 1J The methods described herein, which allow, in part, for the measurement 
of bone structure are useful in both the diagnosis and treatment of osteoporosis. 

5 Ultimately, these techniques could help screen large numbers of women at risk for 
osteoporosis in a highly cost-effective and accurate manner using standard, widely 
available radiographic equipment without the need for expensive dedicated capital 
equipment. It is clear that a program of this type would be powerfully enabling for 
therapeutic intervention with new anabolic or anti-resorptive drugs that are needed to 

1 0 prevent the expected pandemic of osteoporotic fractu res. 

4.1. Kits 

[0212] The invention also provides kits for obtaining information from images, for 
example for obtaining information regarding bone structure, micro-architecture, macro- 
anatomical and/or biomechanical parameters from an image such as a radiograph. In 

15 certain embodiments, the kit comprises one or more computer (e.g., software) 

programs, for example for receiving, analyzing and generating reports based on the 
image(s). In further embodiments, the kits can include calibration phantoms, for 
example calibration phantoms integrated or attachable-to a holder, hygienic cover, x-ray 
film and/or x-ray film holders. 

20 [0213] The invention also provides for therapeutic kits, for example for treating 
osteoporosis or dental disease. In certain embodiments, the kits comprise a calibration 
phantom for use with one or more x-ray films, a computer software product, a database, 
a therapeutic drug and, optionally, instructions for use (e.g., instructions regarding 
positioning the calibration phantom while taking the x-ray, using the software to analyze 

25 the x-ray, dosages and the like. The therapeutic drug can be, for example, anti- 
resorptive or anabolic. 

4.2. Diagnosis and Prediction 
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[0214] In yet another aspect, methods of diagnosing or predicting bone-related 
disorders (e.g., osteoporosis, Paget's Disease, osteogenesis imperfecta, bone cancers), 
periodontal disease or oral implant failure in a subject are provided, for example using 
any of the kits, methods and/or devices described herein. It will be apparent that these 
5 methods are applicable to any bone-relsted disorder including, for example, 

osteoporosis, bone cancer, and the like, as well as to periodontal disease and implant 
failure. 

[02 ia'j Osteoporosis alone is a major public health threat for 25 million 
postmenopausal women and 7 million men. In 1995, national direct expenditures for 

10 osteoporosis and related fractures were $13 billion. Changing demographics, with the 
growth of the elderly population, steadily contribute to increasing numbers of 
osteoporotic fractures and an incipient and potentially economically unmanageable 
epidemic of osteoporosis. Projections put the total cost of osteoporosis in the United 
States alone at more than 240 billion dollars per year in 40 years. 

1 5 [021 6] Less than 20% of the patients know they have the disease and many 
fewer receive physician directed specific therapy. A major impediment in successfully 
dealing with the impending osteoporosis epidemic is not a lack of treatment modalities 
but the inability to identify persons at risk and who require treatment. The limited 
access to osteoporosis testing is largely the result of the high cost of the currently 

20 available systems resulting in a small installed base limited to hospitals and specialty 
clinics. 

[0217] The devices and methods described herein address these and other 
issues by providing inexpensive and reliable bone structural analysis screens and 
resulting diagnosis of bone condition and/or presence of disease. Indeed, while 
25 measurements of bone mineral density (BMD) are technically relatively easy to perform, 
low BMD accounts for considerably less than 100% of fracture risk although it is well 
established that progressive disruption of trabecular structure and architecture 
contribute in a major way to fracture risk in older individuals. 
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[0218] Thus, in certain embodiments, the methods comprise using a computer 
program to analyze bone mineral density or bone structure and/or macro-anatomical 
and/or biomechanical parameters of an image (e.g., x-ray image) and comparing the 
value or measurement obtained from the image with a reference standard or curve, 

5 thereby determining if the subject has a bone-related condition such as osteoporosis or 
thereby determining a subject's fracture risk. The image can also include a calibration 
phantom, for example a calibration phantom as described herein. 
[0219] In certain embodiments, measurements of bone structure can be 
combined or correlated with measurements of macro-anatomical and/or biomechanical 

10 parameters (e.g., cortical thickness on a hip x-ray), for example using statistical or 
mathematical methods, to create an index for the severity of the disease. 
Subsequently, the index can be used for diagnosing osteoporosis or for predicting 
fracture risk combining at least two or more of these bone structure or morphological 
parameters. 

15 4.3. Treatment 

[0220] The methods and devices described herein can also be used to develop 
an appropriate treatment regime for a subject in need thereof. Additionally, the 
invention allows for the ongoing analysis of the efficacy of a subject's treatment regime. 
[0221] Although estrogen deficiency after menopause is one of the most well 

20 documented causes of osteoporosis that can be prevented by hormone replacement 
therapy (HRT), HRT may also cause an increase (approximately 35%) in the risk of 
breast cancer in long-term users. Lancer (1997)350:1047-1059. Consequently, much 
effort has been devoted to developing alternative treatments for osteoporosis. Among 
those treatments, bisphosphonates are becoming increasingly recognized as the 

25 treatment of choice. Lin (1996) Bone 18:75-85; Liberman et al. (1995) N Engl J Med 
333:1437-1443; Mortensen et al. (1998) J Clin Endocrinol Metab 83:396-402. Another 
new class of therapeutic agents recently introduced is the selective estrogen receptor 
modulators (SERMs). Delmas et al. (1997) N Engl J Med 337: 1 64 1 -1 647; Lufkin et al. 
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(1998) J Bone Min Res 13:1747-1754. Anabolic therapies such as parathyroid hormone 
have also been suggested for treatment of osteoporosis. Roe et al. (1999) J Bone Miner 
Res 14(suppl1):S137, Abst#1019; Lane et al. (1998) J Clin Invest 102:1627-33. 
[0222] The combined results of these and other studies suggest that effective 

5 treatments for osteoporosis can be devslopsd once the condition is diagnosed. For 
instance, using any of the methods, kits, and/or devices described herein, the presence 
of osteoporosis in a subject can be diagnosed and that subject provided with 
appropriate therapy (e.g., one or more anti-resorptive agents and/or one or more 
anabolic agents). Periodontal disease can be similarly diagnosed and treatments 

10 ranging from oral hygiene practices to surgery can be recommended. Over time, the 
methods and compositions described herein can be used to assess the efficacy of the 
selected treatment and the treatment regime altered as necessary. For example, a 
subject can be given a one-time or ongoing therapy and images evaluated after such 
therapy to monitor its effectiveness. Thus, in certain embodiments, treatment or 

15 monitoring of treatment of bone related disorders are provided. 
4.4. Decision Trees 

[0223] Thus, diagnosing, predicting, developing treatment regimes, assessing 
treatment efficacy and the like can be readily accomplished using the methods 
described herein. In certain aspects, these applications will be accomplished using 

20 algorithms or decision trees (also known as logic trees or flow charts). One exemplary 
decision tree is provided in regard to predicting bone problems. It will be readily 
apparent that such decision trees are equally applicable to other applications (e.g., 
designing treatment regimes, assessing treatment efficacy, etc.). 
[0224] One exemplary method for predicting bone problems (e.g., osteoporoses, 

25 etc.), periodontal disease or oral implant failure employs a decision tree (also called 
classification tree) which utilizes a hierarchical evaluation of thresholds (see, for 
example, J.J. Oliver, et. al, in Proceedings of the 5th Australian Joint Conference on 
Artificial Intelligence, pages 361-367, A. Adams and L. Sterling, editors, World Scientific, 
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Singapore, 1992; D.J. Hand, etal., Pattern Recognition, 31(5):641-650, 1998; J.J. 
Oliver and D.J. Hand, Journal of Classification, 13:281-297, 1996; W. Buntine, Statistics 
and Computing, 2:63-73, 1992; L. Breiman, etal., "Classification and Regression Trees" 
Wadsworth, Belmont, CA, 1984; C4.5: Programs for Machine Learning, J. Ross 

5 Quinlan, The Morgan Kaufmann Series in Machine Learning, Pat Langley, Series 
Editor, October 1992, ISBN 1 -55860-238-0). Commercial software for structuring and 
execution of decision trees is available (e.g., CART (5), Salford Systems, San Diego, 
CA; C4.5 (6), RuleQuest Research Pfy Ltd., St Ives NSW Australia) and may be used in 
the methods of the present invention in view of the teachings of the present 

10 specification. A simple version of such a decision tree is to choose a threshold bone 
structure and/or macro-anatomical and/or biomechanical or bone mineral density 
reading at a particular anatomical landmark (e.g., edge of mandible or maxilla, the end 
of a tooth root, etc.). If a value is equal to or below the threshold bone data value, then 
more of the image is evaluated. If more of the image is below the threshold value, then 

15 a bone problem, periodontal disease or implant failure is predicted. 

[0225] For example, a first level decision is made by the algorithm based on the 
most recent x-ray images obtained and analyzed as described herein is compared to 
initial thresholds that may indicate an impending or current bone- or periodontal-related 
event. For example, the algorithm may compare the current bone structure 

20 measurements (time=n) or a predicted bone structure measurement (time=n+1) to a 
threshold value. If the bone structure measurement is greater than the threshold value 
then a decision is made by the algorithm to suggest further future x-rays. If the bone 
structure measurement is less than or equal to the threshold level(s) then the algorithm 
continues with the next level of the decision tree. 

25 [0226] The next level of the decision tree may be an evaluation of the subject's 
age and/or gender at time (n) that x-ray is taken, which is compared to a threshold bone 
measurement for "normal" subjects of that age and/or gender. For example, if the 
subject's bone measurement is greater than the threshold bone structure level for that 
particular age and/or gender, then a decision is made by the algorithm to prompt further 
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monitoring in the future. If the information on bone structure is less than or equal to the 
threshold, then the algorithm continues with the next level of the decision tree. 
[0227] The next level of the decision tree may be, for example, an evaluation of 
the subject's soft tissue {e.g., gum) thickness (n), which is compered to a threshold 
5 measurement. For example, if the soft tissue is significantly below or above the normal 
range of thickness, then a decision is made by the algorithm to examine more of the x- 
ray image or to predict a bone-related problem. 

[0228] The decision tree could be further elaborated by adding further levels. For 
example, after a determination that a bone and/or periodontal events are possible, the 
10 subject can be x-rayed again to see if values have changed. Again, age, gender, 

weight, soft tissue thickness and the like can also be tested and considered to confirm 
the prediction. 

[0229] In such decision trees, the most important attribute is typically placed at 
the root of the decision tree. In one embodiment of the present invention the root 
15 attribute is the current bone structure measurement(s). In another embodiment, a 
predicted bone structure measurement at a future time point may be the root attribute. 
Alternatively, bone mineral density and/or implant structure could be used as the root 
attribute. 

[0230] Further, thresholds need not (but can) be established a priori. The 
20 algorithm can learn from a database record of an individual subject's readings and 

measurements. The algorithm can train itself to establish threshold values based on the 
data in the database record using, for example, a decision tree algorithm. 
[0231] Further, a decision tree may be more complicated than the simple 
scenario described above. For example, if soft tissue of a particular subject is very 
25 thick, the algorithm may set a threshold for the bone measurements that is higher or 
lower than normal. 
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[0232] By selecting parameters (e.g., current or future bone information, etc.) and 
allowing the algorithm to train itself based on a database record of these parameters for 
an individual subject, the algorithm can evaluate each parameter as independent or 
combined predictors of disease and/or implsnt failure. Thus, the prediction model is 

5 being trained and the algorithm determines what parameters are the most important 
indicators. A decision tree may be learnt in an automated way from data using an 
algorithm such as a recursive partitioning algorithm. The recursive partitioning 
algorithm grows a tree by starting with all the training examples in the root node. The 
root node may be "split," for example, using a three-step process as follows. (1) The 

10 root node may be split on all the attributes available, at all the thresholds available (e.g., 
in a training database). To each considered split a criteria is applied (such as, GINI 
index, entropy of the data, or message length of the data). (2) An attribute (A) and a 
threshold (T) are selected which optimize the criteria. This results in a decision tree 
with one split node and two leaves. (3) Each example in the training database is 

15 associated with one of these two leaves (based on the measurements of the training 
example). Each leaf node is then recursively split using the three-step process. 
Splitting is continued until a stopping criteria is applied. An example of a stopping 
criteria is if a node has less than 50 examples from the training database that are 
associated with it. 

20 [0233] In a further embodiment, at each level of the decision in the decision tree, 
the algorithm software can associate a probability with the decision. The probabilities at 
each level of decision can be evaluated (e.g., summed) and the cumulative probability 
can be used to determine whether disease and/or implant failure is predicted. 

Receiver Operating Characteristic (ROC) curve analysis can be applied to 

25 decision tree analysis described above. ROC analysis is another threshold optimization 
means. It provides a way to determine the optimal true positive fraction, while 
minimizing the false positive fraction. A ROC analysis can be used to compare two 
classification schemes, and determine which scheme is a better overall predictor of the 
selected event (e.g., evidence of osteoporosis); for example, a ROC analysis can be 
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used to compare a simple threshold classifier with a decision tree. ROC software 
packages typically include procedures for the following: correlated, continuously 
distributed as well as inherently categorical rating scale data; statistical comparison 
between two binomial ROC curves; maximum likelihood estimation of binormal ROC 

5 curves from set of continuous as well as categorical data; and analysis of statistical 
power for comparison of ROC curves. Commercial software for structuring and 
execution of ROC is available (e.g., Analyse-lt for Microsoft Excel, Analyse-lt Software, 
Ltd., Leeds LS12 5XA, England, UK; MedCalc®, MedCalc Software, Mariakerke, 
Belgium; AccuROC, Accumetric Corporation, Montreal, Quebec, CA). 

10 [0234] Related techniques that can be applied to the above analyses include, but 
are not limited to, Decision Graphs, Decision Rules (also called Rules Induction), 
Discriminant Analysis (including Stepwise Discriminant Analysis), Logistic Regression, 
Nearest Neighbor Classification, Neural Networks, and Naive Bayes Classifier. 
[0235] All of these aspects of the invention can be practiced separately or in 

15 combination. Typically, the use of combinations of the embodiments listed above is 
more advantageous. Further, although preferred embodiments of the subject invention 
have been described in some detail, it is understood that obvious variations can be 
made without departing from the spirit and the scope of the invention. 
EXPERIMENTAL 

20 [0236] Below are examples of specific embodiments for carrying out the present 
invention. The examples are offered for illustrative purposes only, and are not intended 
to limit the scope of the present invention in any way. 
Example 1: In vivo reproducibility and in vivo diagnostic sensitivity 

A. Dental X-Rays 

25 [0237] In order to test in vivo reproducibility of data obtained from dental x-rays, 
the following experiment was performed. Subjects sat in a dental chair and an x-ray 
was taken of the area of the incisor teeth and of the molar teeth of the mandible. A 
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calibration phantom step wedge was attached to the dental x-ray film. The dental x-ray 
film was exposed using standard x-ray imaging techniques for x-rays of the incisor area. 
The subjects walked around for 1 5 minutes at which point that test was repeated using 
the same procedure. 
5 [023S] X-ray films were digitized on a commercial flat-bed scanner with 

transparency option (Acer ScanPremio ST). The regions of interest (ROIs) were placed 
manually at the same position with respect to the dental roots in all digitized x-rays of 
the same subject using the NIH Image software program (http://rsb.info.nih.gov/nih- 
image/Default.html). The reproducibility of the measurement of the average gray values 
10 inside the ROIs was determined as the coefficient of variation (COV=standard deviation 
of measurements/mean of measurements). Overall results are given as root mean 
square (RMS = -^^xfjn) over both subjects. The data are summarized in Table 2. 



Table 2: Reproducibility of measurements of average gray values 



Region 


COV Subject A 


COV Subject B 


RMS 


Incisor 


2.9% (n=3) 


5.9% (n=3) 


4.6% 


Molar 


3.0% (n=3) 


4.1% (n=4) 


3.6% 




All regions: 


4.2% 



[0239] The data show that reproducibility is achieved that is already comparable 
with that of many ultrasound systems to diagnose osteoporosis. 

B. Hip Radiographs 

20 [0240] To test whether bone texture analysis in hip x-rays can detect differences 
between normal and osteoporotic bone, sample hip x-ray images were acquired in two 
patients with a Fuji FOR 5000 computed radiography system (Fuji Medical Systems, 
Stemford, CT). The first patient had normal bone mineral density in the hip as measured 
by DXA. In the second patient, femoral neck BMD measured by DXA was one standard 

25 deviation below normal. 
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[0241] For x-ray imaging, patients were positioned on the x-ray table in supine 
position, parallel to the long axis of the table. The patient's arms were placed alongside 
their body. Patient comfort was ensured with a pillow underneath the patient's neck. 
However, no pillows were used underneath the knees. The x-ray technologist checked 
5 that the patient lies straight on the table by looking from the head down towards the feet 
(which were placed in neutral position with the toes pointing up. The ray was centered 
onto the hip joint medial and superior to the greater trochanter. 

[0242] Anteroposterior hip radiographs were acquired using the following 
parameters: Film-focus distance: 100 cm; tube voltage: 65 kVp; exposure: phototimer 

10 for automatic exposure or approximately 20 mAs for manual exposure; collimation: 
limited to the hip joint, including proximal femoral diaphysis; centering: over femoral 
head (see above); tube angulation: zero degrees. An aluminum step wedge (BioQuest, 
Tempe, AZ) was included in the images to calibrate gray values before further image 
analysis. Processing was performed using ImageJ, a Java version of NIH image 

15 (http://rsb.info.nih.gov/ij/). 

[0243] Six regions of interest were selected manually at the approximate 
locations as shown in Figure 9. Trabeculae were extracted through background 
subtraction. The resulting binarized images are shown in the Figures. In a next step, 
the trabecular bone in the selected regions of interest was skeletonized. 

20 [0244] The binarized ROI's in the normal and the osteopenic patient were used to 
determine the trabecular density ratio (trabecular area vs. ROI area). The following 
bone structure measurements were obtained from the skeletonized ROI's; mean 
segment length, total skeleton length (normalized by ROI area), skeleton segment count 
(normalized by ROI area), and skeleton node count (normalized by ROI area). Results 

25 are shown in Tables 3 through 7. 



Table 3; Trabecular Density Ratio (Trabecular Area / ROI Area) 





ROI A 


ROIB 


ROIC 


ROI D 


ROI E 


ROI F 


Normal 


0.473 


0.482 


0.514 


0.494 


0.476 


0.485 
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ROI A 


ROI B 


ROIC 


ROI D 


ROI E 


ROIF 


Osteopenia 


0.382 


0.455 


0.492 


0.426 


0.424 


0.455 


% Osteopenia vs. Normal 


81% 


94% 


96% 


86% 


89% 


94% 


Table 4: Mean skeleton segment length 




ROI A 


ROI B 


ROIC 


ROI D 


ROIE 


ROI F 


Normal 


7.116 


8.071 


10.765 


8.175 


8.272 


7.313 


Osteopenia 


7.146 


9.877 


10.004 


6.699 


8.607 


9.750 


% Osteopenia vs. Normal 


100% 


122% 


93% 


82% 


104% 


133% 


Table 5; Total Skeleton Length (normalized by ROI area) 




ROI A 


ROIB 


ROIC 


RO I D 


ROI E 


ROI F 


Normal 


0.0736 


0.0758 


0.0906 


0.0889 


0.0806 


0.0785 


Osteopenia 


0.0503 


0.0589 


0.0672 


0.0584 


0.0681 


0.0543 


% Osteopenia vs. Normal 


68% 


78% 


74% 


66% 


84% 


69% 


Table 6: Skeleton segment count (normalized by ROI area) 




ROI A 


ROIB 


ROI C 


ROID 


ROI E 


ROIF 


Normal 


0.0100 


0.0094 


0.0084 


0.0109 


0.0097 


0.0107 


Osteopenia 


0.0070 


0.0060 


0.0067 


0.0087 


0.0079 


0.0056 


% Osteopenia vs. Normal 


68% 


63% 


80% 


80%> 


81% 


52% 


Table 7: Skeleton node count (normalized by ROI area) 




ROI A 


ROI B 


ROI C 


ROID 


ROI E 


ROI F 


Normal 


0.0198 


0.0210 


0.0229 


0.0244 


0.0156 


0.0240 


Osteopenia 


0.0090 


0.0117 


0.0132 


0.0113 


0.0088 


0.0081 


% Osteopenia vs. Norma 


46% 


56% 


58% 


47% 


56% 


34% 



5 

These results demonstrate that the evaluation of trabecular structure reveals significant 
differences between normal and osteopenic bone and that selective analysis of 
trabeculae oriented in certain directions in the different ROI allows for the assessment 
of structures critical for biomechanical stability of the proximal femur. 
10 C. Spine Radiographs 

[0245] To test whether bone texture analysis in spine x-rays can detect 
differences between normal and osteoporotic bone, sample spine x-ray images will be 
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acquired in more than one patient. In the spine, the bone structure parameters can be 
measured in the L1 , L2, L3 and L4 vertebral bodies unless obscured by superimposed 
ribs, iliac crest or bowel gas. The first patient will provide control data provided the 
patient has normal bone mineral density in the spine. In the second patient and 
5 subsequent patients, spine BMD will be measured. 

[0246] Regions of interest will be selected manually at the approximate locations as 
shown in Figure 22. Trabeculae will be extracted through background subtraction. In a 
next step, the trabecular bone in the selected regions of interest was skeletonized. Fig. 
24 depicts an example of an application of structure extraction and measurement for 
10 therapeutic monitoring using spine x-ray. White outline of extracted structure are show 
in (a) before treatment, and (b) after treatment. 

D. Knee/Tibial Radiographs - Osteoporosis 

[0247] To test whether bone texture analysis in knee and tibial x-rays can detect 
differences between normal and osteoporotic bone, sample x-ray images will be 
15 acquired in more than one patients. The first patient will provide control data provided 
the patient has normal bone mineral density or bone structure in the tibia or femur. In 
the second patient and subsequent patients, joint BMD or bone structure will be 
measured. 

[0248] Regions of interest will be selected manually at the approximate locations 
20 as shown in Figure 23. The ROI can, for example, be the region immediately below the 
tibial plateau subchondral bone. Trabeculae will be extracted through background 
subtraction. In a next step, the trabecular bone in the selected regions of interest is 
skeletonized. 

D. Knee/Tibial Radiographs - Arthritis 
25 [0249] To test whether bone texture analysis in knee and tibial x-rays can detect 
differences between normal patients and patients with arthritis, sample x-ray images will 
be acquired in more than one patients. The first patient will provide control data 
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provided the patient has normal bone mineral density or bone structure in the tibia or 
femur. In the second patient and subsequent patients, joint BMD or structure will be 
measured. 

[0250] Regions of interest will be selected manually at the approximate locations 
5 as shown in Figure 23. The ROI can, for example, be the region immediately below the 
tibial plateau subchondral bone. Trabeculae will be extracted through background 
subtraction. In a next step, the trabecular bone in the selected regions of interest is 
skeletonized. 

Example 2: Image Processing Techniques 

10 [0251] Techniques to analyze structure of trabeculae in different regions of the 
femoral head, neck, and proximal shaft are developed in Matlab (The MathWorks, Inc., 
Natick, MA) on PC's. The following techniques (modules) are developed: algorithms for 
software analysis of density, length, thickness, and orientation of trabeculae in different 
regions of interest (ROI) in the radiograph and a technique for automated placement of 

15 these ROI. 

[0252] Six regions of interest are selected in the proximal femur for bone 
microstructure evaluation. The size and shape of these ROI are designed to capture 
the local changes of trabecular density and structure (see, e.g., Figure 9), and may 
reflect the location of the different compressive and tensile groups of trabeculae. Singh 
20 et al. (1970) J Bone Joint Surg Am. 1970. 52:457-467. Thus, a classification scheme 
based on statistical convergence of multiple parameters that would provide a high 
precision index for predicting hip fractures is developed. 
Example 3: Bone Structure Analysis of Hip Radiographs 

[0253] The trabeculae in the femur is extracted using the background subtraction 
25 method, essentially as described in Geraets et al. (1998) Bone 22:165-173. A copy of 
the image is blurred with a 15x15 Gaussian filter, and the result represents the non- 
uniform background. This background image is subtracted from the original image to 
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obtain an image of trabecular structure. This image is then transformed into binary 
image of trabecular structure by applying a threshold value of 0. An example of the end 
result is shown in Figure 10. 

[0254] In a second step, parameters relevant to the geometry and connectivity of 
5 trabecular structure are measured on the trabecular skeleton or centerline. The 
skeletonization is performed using morphological hit-or-miss thinning for example as 
described in Soille, "Morphological image analysis: principles and application" Springer, 
1998: p. 129-154. The branch points and end points of the skeleton network are 
detected, and the skeleton segments are classified as free-end segments and node-to- 
10 node segments. 

[0255] One or more of the following parameters from the binarized and from the 
skeletonized ROI's are used: trabecular density; ratio of trabecular area to total ROI 
area; trabecular perimeter; star volume (Ikuta et al. (2000) J Bone Miner Res. 18:271- 
277; Vesterby (1990) Bone 11:149-155); trabecular bone pattern factor (Hahn et al. 

15 (1992) Bone 13:327-330); Euclidean distance transform; assessment of trabecular 
orientation using Fourier analysis; and orientation-specific trabecular assessment. 
Further, one or more of the following parameters can be measured in each ROI on the 
network of skeletonized trabeculae as a whole, all skeleton segments, and each type of 
segment: segment count; segment length; angle of segment orientation; and 

20 Interconnectivity Index (Legrand et al. (2000) J. Bone Miner Res. 15:13-19): normalized 
ratio of the number of node-to-node segments to free-end segments. 
[0256] For example, in Euclidean Distance Transform each pixel on the binarized 
trabeculae is assigned a value equal to its Euclidean distance from the structure 
boundary. Thus, thicker trabeculae will have larger distance transform values in the 

25 center, thereby estimating trabecular thickness calculates the mean of the distance 
transform values along the trabecular skeleton (see Figure 11). Further, multiplying this 
value by 2 provides a measurement of trabecular thickness. 
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[0257] Similarly, predominant trabeculae orientation may be evaluated using the 
2D Fast Fourier Transform (FFT). A rectangular region is selected within each ROI and 
multiplied with a 2D Kaiser window before applying the transform (see Figure 12, left). 
The log of the Fourier magnitude is taken to fonn an image representing the frequency 
5 domain of the ROI. The result is then filtered with a 5x5 Gaussian filter to reduce local 
variation. An example image is shown in Figure 12, center. The Fourier image is 
subsequently thresholded at a fixed magnitude level. This binary image is resampled to 
a square image to normalize the length of the vertical and horizontal axes, and the 
direction and length of its major axis are determined (Figure 12, right). The angles will 

10 be measured with respect to the axes of the femoral neck and shaft. The axes are 

determined by fitting lines to the two longest segments of the centerline of the binarized 
femur (see also Figure 14). The ROI's are located such that they include the different 
groups of compressive and tensile trabeculae in the proximal femur that each can be 
characterized by a specific direction. A fully automated technique to evaluate the 

15 different quantitative structural parameters explained above for those trabeculae in each 
of the ROI that are oriented in the characteristic direction expected for the particular ROI 
is developed. 

[0258] The orientation of each trabecular skeleton segment is found through the 
gradient of the line fitted to the skeleton points. Based on this orientation information, 
20 only those trabeculae are considered in the evaluation of the structure parameters that 
are approximately oriented in the characteristic direction for a particular ROI. 

[0259] As will be appreciated by those of skill in the art, all measurements can be 
constrained by one or more desired orientation by measuring only segments within 
specified angle ranges. The statistics of watershed segments include: number of 
25 segments, total area of segments, average area of segments, standard deviation of 

segment area, smallest segment area, and largest segment area. These segments are, 
however, general in nature. 
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[0260] When evaluating the hip, additional parameters can be considered. 
Parameters include, for example, shaft angle, neck angle, diameter of the femur neck, 
the hip axis length, the largest cross-section of the femur head, the average thickness of 
the cortical region within a ROI, the standard deviation of cortical thickness within a 

5 ROI, or the maximum or minimum thickness of the cortical thickness within a ROI. 
[0261] In contrast, when evaluating the spine, additional parameters to be 
considered include, for example, all parameters on vertical structures, all parameters on 
horizontal structures, vertebral cortical thickness, maximum vertebral height, minimum 
vertebral height, average vertebral height, anterior vertebral height, medial vertebral 

10 height, posterior vertebral height, maximum inter-vertebral height, minimum inter- 
vertebral height, and average vertebral height. 

[0262] The knee and tibial region can be evaluated using the additional 
parameters of: average medial joint space width, minimum medial joint space width, 
maximum medial joint space width, average lateral joint space width, minimum lateral 

15 joint space width and maximum lateral joint space width. 

[0263] As will be appreciated by those of skill in the art, the additional parameters 
listed for these exemplar anatomies above can include other parameters. Additionally, 
parameters can be evaluated for other anatomies not specifically set forth without 
departing from the scope of the invention. 

20 Example 4: Multidimensional Classification 

[0264] Example 3 describes a number of parameters that are measured to 
assess trabecular structure in different regions of the proximal femur. In this Example, 
the different structural parameters are combined in each section, and a single index is 
determined over all regions of interest. 

25 [0265] A training set of hip x-ray images of a group of subjects are divided into 
the two categories "osteoporosis" and "no osteoporosis", based on previous DXA 
results. Subsequently, for all x-rays in the training set, the parameters listed in Example 
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3 are calculated for all regions of interest placed as described in Example 3, resulting in 
a set of m-dimensional prototype feature vectors f t = {f n ,-,f im f for the training set 

I = {/,},* = 1,...,R. 

[0266] For each parameter a single scalar index value is calculated. All index 
5 values are combined into one n-dimensional feature vector. In one step, the system is 
trained with the data from clinical validation studies with premenopausal, 
postmenopausal healthy and postmenopausal osteoporotic subjects. The subject 
groups are preferably divided into a "fracture" and a "no fracture" category. The feature 
vectors calculated from the x-ray images are used as prototype patterns. 
10 [0267] For each patient, a feature vector is calculated from the x-ray as 

calculated for the prototype patterns and an individual patient classified as category C if 
the majority of the k closest prototype patterns is of the category C. The distance d 
between the patient's feature vector f = (f l ,f 2 ,-,f„) T and a prototype pattern 
p = (P\>P2>--->Pnf ' iS defined by the Euclidean norm l_ 2 : 

15 d (/, p) - L 2 (/, p) = (/ ; -Pjj 1 

[0268] The optimum scale for the different parameters is also preferably 
determined. However, for some parameters differences in the index values between 
the categories is smaller than for others. Also, the optimum k will be determined. 
Increasing k is expected to improve the accuracy of the classification, but it has to be 

20 smaller than the number of prototypes in each category. The exact percentage value of 
the majority of the k closest prototype patterns that determines the classification 
provides a measure for the reliability of the classification. The higher the percentage of 
prototype patterns from a particular category C, the more significant the information 
provided by the classification is likely to be. 

25 [0269] This classification approach is validated with a series of leave-one-out 
experiments using the 0° neutral position images of the femoral position study (see 
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Example 8) and the baseline hip x-rays of the short-term in vivo reproducibility study. 
For these experiments, each subject is preferably used as a test case once. The 
training set for the system consists of the patterns calculated for all or most of the 
remaining subjects. The test case is correctly classified using this training set, and the 
diagnostic sensitivity and specificity of the combination of bone structure parameters is 
determined. 

[0270] In addition to the measurements described above (which provide index 
values for the parameters "length of trabeculae", "direction of trabeculae and 
anisotropy", and "trabecular thickness"), additional measurements for other parameters 
in the classification system that have been explored in the past to study bone density 
and structure from x-ray, CT, and MR images such as: (1) mean pixel intensity; (2) 
variance of pixel intensity; (3) Fourier spectral analysis; (4) fractal dimension; (5) 
morphological parameters such as the trabecular area, trabecular periphery, total 
trabecular length, number of terminal and branch points, as well as similar parameters 
for the bone marrow can be used. 

Example 5: Automated Placement of Region of Interest (ROI) 
[0271] Analysis of x-rays (e.g., hip radiographs) may be facilitated by 
development of techniques that locate one or more regions of interest (ROI) used for 
the calculation of the structural parameters of the trabecular bone. For example, the 
general position of the femur can be located using a binary image of the hip radiographs 
thresholded at the appropriate gray value. In a typical hip radiograph, the femur is a 
bright structure extending from the pelvis. (Figure 13). By thresholding the digitized 
radiograph at the typical femur intensity value, a binary image showing the femur is 
produced. The relatively thin structure of the femoral shaft can be extracted by applying 
a morphology operation on the binary image. The morphological top-hat filter (opening 
subtracted from input) with an upright rectangular structuring element segments the 
femoral shaft. The result is shown in Figure 13 with outline of the binarized femur 
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superimposed on the original radiograph. The region is cropped for further processing, 
preferably leaving enough room to include the femoral head. 

[0272] To position the set of predetermined ROI, a regularized active shape 
algorithm can be used (Behiels et al. (1999) Proceedings of the 2nd International 
5 Conference on Medical Image Computing and Computer-Assisted Intervention - 
MICCAI'99, Lecture notes in Computer Science 1679:128-137; Cootes (1994) Image 
and Vision Computing 12:355-366). A general model of the proximal femur is created by 
manually outlining the shape in a training set of typical hip radiographs to form a mean 
shape. The six predefined ROI are then embedded into this model. This mean model is 

10 scaled down 80%, isometrically along its centerline. This transformation is applied to the 
predefined ROI as well. The outline of the rescaled model is then used as the initial 
template and is positioned within the proximal femur in the input image. The control 
points of the contour are subsequently expanded outwards away from the nearest 
centerline point. The energy function to be optimized in this iterative process can take 

15 into account local features, such as gradient, intensity, deviation from the mean model, 
and curvature of contour segments. Figure 14 illustrates the propagation of the initial 
control points towards the femur edge. When the iteration is completed, a deformation 
field for the model area is calculated. This deformation field is interpolated for the 
model ROI inside the boundaries of the femur model. The result is a new set of ROI that 

20 is adapted to the input image, but similar to the model ROI with respect to anatomical 
landmarks (see Figure 9). 

Example 6: Data Analysis 

[0273] Patients are selected into one of three groups: healthy premenopausal 
(PRE); healthy postmenopausal (POST), and osteoporotic postmenopausal (OSTEO) 
25 women. All groups are studied by: (1 ) dental x-ray images of the periapical and canine 
region; (2) quantitative computed tomography of the spine and (3) hip; (4) dual x-ray 
absorptiometry of the spine and (5) hip; (6) single x-ray absorptiometry of the 
calcaneus, and (7) ultrasound of the calcaneus using standard techniques. A diagnosis 
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of osteoporosis is made when at least one atraumatic vertebral fracture as determined 
by a semi-quantitative assessment of morphologic changes of the thoracic and lumbar 
spine on lateral conventional x-rays is observed. 

[0274] The means and standard deviations of the different bone structure 
5 measurements (see above) and bone mineral density measurements (mandibular BMD, 
QCT spine, QCT hip, DXA spine, DXA hip, SXA calcaneus, ultrasound calcaneus) are 
calculated for each patient group. The Student's t-test (t-values and p-values) and 
percent decrement are used for comparing the different measurements for reflecting 
intergroup differences. Annual, age-related changes are expressed as percent changes 

10 relative to the predicted values at age 30 and as fractional standard deviation (SD) of 
PRE. Correlations with age along with p-values are also be reported. Odds ratios (for 
1SD change in the measured parameter) and 95% confidence limits based on the age- 
adjusted logistic regression are calculated to measure the discriminative ability (for 
discriminating between the postmenopausal osteoporotic and the normal 

15 postmenopausal group) and the risk of osteoporotic fracture associated with the 
measured parameter. The pairwise comparisons of the discriminative abilities are 
tested using age-adjusted receiver operating characteristic (ROC) curve analysis. 

[0275] Pairwise comparisons of all techniques are obtained by pooling all 
subjects (PRE, POST, OSTEO) and using Pearson's correlation coefficients (r), percent 
20 standard errors of the estimate (CV), and p-values for testing significance of 
correlations. 

[0276] To compare measurements for their diagnostic ability, a kappa score 
analysis is performed on the normal postmenopausal women (POST) and the 
osteoporotic postmenopausal women (OSTEO). This is done by classifying every 
25 woman from the postmenopausal groups as osteopenic if her T-score with respect to 
the reference group (PRE) is less (or in case of structural parameters also greater) than 
2.5. The T-score for an individual woman and a particular measurement is defined as 
the measurement minus the mean measurement of young normals (PRE) divided by the 
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SD of the measurement in the PRE group. Note that the T-score is measuring the 
position of an individual woman with respect to the PRE group and is different from the 
Student's t-value. 

Example 7: Longitudinal Monitoring of Bone Structure 

5 [0277] Algorithms and software to match follow-up dental x-rays obtained at a 
time point T 2 relative to baseline x-rays of the mandible obtained at an earlier time point 
Ti are developed. For purposes of monitoring of therapeutic response, bone structure 
parameters have to be measured at the same location of the mandible at different 
points in time. Thus, in order to compensate for differences in patient positioning and in 
10 order to find corresponding regions of interest (ROI's) for comparison of the results 

between baseline and follow-up examinations, it is desirable to register two dental x-ray 
images. 

[0278] Due to possible slight differences in the projection angle of the x-ray beam 
on the film in the two images to be registered, an elastic matching step is preferably 
15 included. The first step, however, is a global affine transformation, for which the mutual 
information is used as a cost function. Wells et al. (1996) Medical Image Analysis 1:35- 
51. The mutual information I m ,n of two images M and N is defined as 

[0279] Here, the gray values occurring in the two images are regarded as random 
20 variables, and the mutual information provides a measure of the strength of the 
dependence between these variables. p M and p N are the distributions of M and N 
respectively, and p M N is the joint distribution of M and N. Maintz et al. (1998) SPIE 
Medical Imaging - Image Processing. These distributions can be approximated from the 
marginal and joint gray value histograms, more accurately with the use of a Parzen 
25 window function. Powell's method can be used as an optimization scheme to find the 
best affine transformation for N to match it with M. Press et al. ("Numerical Recipes in 
C." 2nd edition, 1992, Cambridge University Press. 
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[0280] This global transformation is followed by local elastic adjustments to 
improve the match. To achieve this, the conditional probability densities p(n|m) are 
estimated from the joint histogram of the globally registered images. The transformation 
vector field t(x) is then determined such that N(x-i(x)) is as similar to m(\) as possible 
5 by maximizing the local gray value correspondence, which for a fixed value of x is 
defined as 

c x (t) = jVx' - x)p(N(x' - 1) | M(x'))dx' . 

[0281] Here, w is a window function whose width determines the size of the 
region that is used to compute t(x). To determine the window function, an approach 
10 similar to the one described in Warfield et al. "Brain Warping" 1 999, Academic Press, 
p:67-84 is used. A number of successively wider window functions Wj are combined into 
a single window w='^W i w l , where the weights W, are given as 

W t ==r-i det(e,) with Q, = fw j (^'-x)V7V(x')ViV 7 '(x')^' . 

[0282] The exact location of the ROI after automatic placement in the baseline 
15 image for a particular patient is kept in a database. When the patient returns for a 
follow-up exam, the new image is registered with the baseline image, and thus 
transformed into the coordinate system of the baseline image. The bone structure in the 
registered follow-up x-ray can then be measured at exactly the same position as in the 
baseline image. 

20 Example 8: Influence of Positioning of the Femur on Bone Structural 
Measurements 

[0283] The effect(s) of the positioning of the femur on each parameter of the 
bone structure assessments is (are) examined. Hip x-rays are obtained in normal 
postmenopausal women and postmenopausal women with osteoporosis in neutral 
25 position and in various degrees of internal and external rotation. 
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[0284] The diagnosis of osteoporosis is made when at least one atraumatic 
vertebral fracture as determined by a semi-quantitative assessment of morphologic 
changes of the thoracic and lumbar spine on lateral conventional radiographs is 
observed. See, also, Genani et al. (19Q3) J. Bone Miner Res. 8:1 137-1 148. 

5 [0285] Standard anteroposterior hip radiographs are obtained with the extremity 
at 30° internal rotation, 15° internal rotation, 0°, 15° external rotation, and 30° external 
rotation. These angles are achieved by placing the foot and ankle against a 30° or a 15° 
degree wedge in either internal or external rotation of the femur. The foot is secured 
against the wedge using Velcro straps. 

10 [0286] The effect of positioning is assessed by calculating the pair wise 
coefficient of variation (CV%) between the results for the 0° position and the other 
positions for each individual subjects. The angular dependency will be expressed for 
each of the angles 30° internal rotation, 15° internal rotation, 15° external rotation, and 
30° external rotation as the root-mean-square of these CV% values over all subjects. In 

15 general, parameters with the least dependency on angular positioning of the femur are 
selected. 

[0287] If the pair wise coefficient of variation between the results for the 0° 
neutral position and the 15° internal or external rotation position exceed 10% for the 
majority of the structural parameters measured, a foot holder that fixes the patients' foot 

20 in neutral position can be used The foot holder is designed with a base plate extending 
from the mid to distal thigh to the heel. The base plate preferably sits on the x-ray table. 
The patients' foot is positioned so that the posterior aspect of the heel is located on top 
of the base plate. The medial aspect of the foot is placed against a medial guide 
connected rigidly to the base plate at a 90° angle. A second, lateral guide attached to 

25 the base plate at a 90° angle with a sliding mechanism will then be moved toward the 
lateral aspect of the foot and will be locked in position as soon as it touches the lateral 
aspect of the foot. The foot will be secured to the medial and lateral guide using Velcro 
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straps. It is expected that the degree of involuntary internal or external rotation can be 
limited to less than 5° using this approach. 

Example S: Influence of X-Ray Tube Angulation on Bone Structural 
ries'furemei-ite 

5 [0288] The effect(s) of the positioning of the x-ray tube on each parameter of the 
bone structure assessments is (are) examined. Dental x-rays are obtained in normal 
postmenopausal women and postmenopausal women with osteoporosis. The diagnosis 
of osteoporosis is made when at least one atraumatic vertebral fracture as determined 
by a semi-quantitative assessment of morphologic changes of the thoracic and lumbar 
10 spine on lateral conventional radiographs is observed. See, also, Genant et al. (1993) 
J. Bone Miner Res. 8:1137-1148. 

[0289] Standard anteroposterior dental radiographs are obtained in the incisor 
region of the mandible. The x-ray tube is aligned with an angle of 0°, 10°, 20°, 30°, and - 
10°, - 20°, and -30° relative to the dental x-ray film. These angles are achieved with use 
15 of a goniometer applied to the metal tube located in front of the dental x-ray tube. The 
dental x-ray film is positioned at the posterior mandibular wall in the incisor region. 

[0290] The effect of positioning is assessed by calculating the pair wise 
coefficient of variation (CV%) between the results for the 0° position and the other tube 
positions for each individual subject. The angular dependency will be expressed for 
20 each of the angles as the root-mean-square of these CV% values over all subjects. 

[0291] The results indicate that a 10 degree tube angulation can result in a 12% 
error in apparent density. 

[0292] A mechanical alignment system is then applied to the Rinn holder. For 
this purpose, an extension tubing is attached to the Rinn holder. The extension tubing 
25 is designed so that its inner diameter is slightly greater (and fits over) than the outer 
diameter of the dental x-ray system metal tube (Fig. 15). The dental x-ray system metal 
tube is then inserted into the extension tubing attached to the Rinn holder that reduces 
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alignment error of the x-ray tube relative to the x-ray film. One group of patients then 
undergo two x-rays each of the incisor region. The results indicate that the short-term 
in-vivo reproducibility error of dental bone density and bone structure measurements is 
reduced with use of the mechanical alignment system by reducing x-ray tube angulation 
5 relative to the dental film and the anatomic landmarks in the mandible. 

Example 10: Measurement of Bone Density, Bone Structure, Macro-Anatomical 
Parameters and Biomechanical Parameters and Selecting Therapy 
[0293] An x-ray image of a mandible or a hip or spine or other bone is analyzed 
using a computer program capable of assessing bone density, bone structure, macro- 

10 anatomical parameters, or biomechanical parameters, for example as described above. 
The computer program derives a measurement of one or more bone density, bone 
structure, macro-anatomical or biomechanical parameters of the trabecular bone. The 
measurement of the parameter(s) is compared against a database containing 
information on said one or more parameters in normal, healthy age-, sex-, and race 

15 matched controls. If the patient's measurement differs by more than 2 standard 

deviations from the age-, sex-, and race matched mean of normal, healthy subjects, a 
report is sent to the physician who then selects a therapy based on the 
measurement(s). 

Example 11: Measurement of Bone Density, Bone Structure, Macro-Anatomical 
20 Parameters and Biomechanical Parameters and Monitoring Therapy 

[0294] One or more x-ray images (mandible, hip or spine or other bone) are 
obtained from a patient undergoing therapy for osteoporosis, for example using an 
anabolic or an antiresorptive drug at two different time points T1 and T2. The x-rays are 
analyzed using a computer program capable of assessing bone density, bone structure, 
25 macro-anatomical parameters, or biomechanical parameters. The computer program 
derives a measurement of one or more parameters of the bone for both time points T1 
and T2. The measurement of the bone density, bone structure, macro-anatomical, or 
biomechanical parameter(s) at T1 and T2 is compared against a database containing 
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information on said one or more parameters in normal, healthy age-, sex-, and race 
matched controls for each time point. If the results indicate that the patient has lost 5% 
or more bone between time points T1 and T2 despite therapy, a physician selects a 
different, more aggressive therapy. 
5 Example 12: Measurement of i"iaci'c-i n; comic -. I and/or biomechanics! Parameters 
[0295] A hip radiograph is obtained using standard techniques and including a 
calibration phantom as described herein. The reference orientation of hip x-rays is the 
average orientation of the femoral shaft. 

A. Edge-detection 

10 [0296] A global gray level thresholding is performed using a bi-modal histogram 
segmentation algorithm on the hip x-ray generates a binary image proximal femur. 
Edge-detection of the hip x-ray can be used. Optionally, edge-detection methods are 
further refined by obtaining breaking edges detected into small segments and 
characterizing the orientation of each segment, thereby obtaining the outline of proximal 

15 femur. Each edge segment is then referenced to a map of expected proximal femur 
edge orientation and to a map of probability of edge location. Edge segments that do 
not conform to the expected orientation or are in low probability regions are removed. 
Morphology operations are applied onto the edge image to connect edge 
discontinuities. The edge image forms an enclosed boundary of the proximal femur. 

20 The region within the boundary is then combined with the binary image from global 
thresholding to form the final mask of the proximal femur. 

[0297] Within a selected region of interest, edge detection is applied. 
Morphology operations are applied to connect edge discontinuities. Segments are 
formed within enclosed edges. The area and major axis length of each segments are 
25 then measured. The regions are also superimposed on the original gray level image 
and the average gray level within each region is measured. The cortex is identified as 
the segments that are connected to the boundary of the proximal femur mask, that has 



WO 2004/086972 



PCT7US2004/00<M<>5 



the greatest area, longest major axis length and has a mean gray level above the 
average gray level of all enclosed segments within the proximal femur mask. 

The segment identified as cortex is then skeletonized. The orientation of the Gortex 
skeleton is verified to conform to the orientation map of proximal femur edge. Euclidian 
distance transform is applied to the binary image of the segment. The values of distance 
transform value along the skeleton are sampled and statistics (average, standard 
deviation, minimum, maximum and mod) measured. 

[0298] As will be appreciated by those of skill in the art, measurements of macro- 
anatomical parameters described here can be applied to hip, spine or knee radiographs 
with modifications to adapt to the shape, scale and location of macro-anatomical 
features specific to the anatomical region. 
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CLAIMS 

1 . A method to derive information regarding one or more bone parameters from 
an image comprising: 

(a) obtaining an image comprising bone from a subject; 
5 (b) defining two or more regions of interest (ROIs) in the image; and 

(c) analyzing a plurality of positions in the ROIs to determine one or more 
parameters selected from the group consisting of bone microarchitecture, bone macro- 
anatomy, biomechanics! parameters and combinations thereof of the ROIs. 

10 2. The method of claim 1, wherein the ROIs are overlapping. 

3. The method of claim 1, wherein the positions analyzed in the ROIs are at 
regular intervals in the image. 

15 4. The method of claim 1 , wherein the positions analyzed in the ROIs are at 

irregular intervals in the image. 

5. The method of claim 1, wherein the parameter is bone micro-architecture and 
the positions analyzed are at regular intervals. 

20 

6. The method of claim 1 , wherein the parameter is bone macro-anatomy and 
the positions are analyzed are at irregular intervals in the image. 

7. The method of any of claims 1-6, wherein the image is two-dimensional. 

25 

8. The method of claim 7, wherein image is an x-ray image. 

9. The method of any of claims 1-6, wherein the image is three-dimensional. 
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10. The method of any of the preceding claims wherein the image is an 
electronic image. 



11. The method of claim 1, wherein the subject is an osteoporosis subject. 



5 



12. A method of generating a map of one or more bone parameters, comprising 

(a) obtaining information on bone parameters according to the method of any of 
claims 1-11; and 

(b) identifying regions of the image that exhibit similar parameter characteristics, 
10 thereby creating a parameter map of the image. 

13. A method of predicting a fracture path in a subject, comprising: 

(a) generating multiple parameter maps according to the method of claim 12; 

(b) generating a composite parameter map from the multiple parameters maps of 
15 step (a); and 

(c) analyzing the composite parameter map to identify possible fracture paths. 

14. A method of predicting a fracture path in a subject comprising: 

(a) analyzing one or more parameter maps preparing according to the method of 
20 claim 12, wherein the analysis is watershed segmentation analysis or Markov random 
field analysis', and 

(c) identifying possible fracture paths based on the analysis of step (a), thereby 
predicting a fracture path in the subject. 

25 1 5. A method of predicting the risk of fracture in a subject comprising: 

(a) generating a finite element model from one or more parameter maps obtained 
according to the method of claim 12; 

(b) applying simulated force vectors that would occur during a fracture incident to 
the finite element model generated in step(s); and 
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(c) determining the minimum forces required for fracture to occur, thereby 
estimating the risk of fracture. 

18. A method of determining the risk of fracture in a subject comprising: 

(a) predicting a fracture path according to the method of claim 13 or claim 14; 

(b) evaluating one or more selected bone parameters along the predicted 
fracture path, thereby estimating the risk of fracture. 

17. A method of treating a subject with bone disease comprising 

(a) obtaining an image from a subject; 

(b) analyzing the image obtained in step (a) as described in any of claims 1-11; 

(c) diagnosing a bone disease based on the analysis of step (b); and 

(d) selecting and administering a suitable treatment to said subject based on said 
diagnosis. 
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FIG„ 12 



FIG. 13 
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Figure 16. An example of a regular interval sampling field for microarchitecture (+) 
and a higher density sampling field for macro-anatomical features(*) on a femur 
radiograph. White rectangles are examples of overlapping window positioning. 




Figure 17. Watershed segmentation boundaries superimposed on a parameter map. 
The two white lines are the actual fracture paths resulted from an in- vitro 
mechanical loading test. 
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and window 






• Regular or irregular density sampling grid 
optimized for resolution and efficiency. 

• Window size and shape optimized for 
specificity. 












Overlapping window 
analysis 






• For each sampling point, analyze density, 



microarchitecture and macro-anatomical 
parameters. 

Generate matrices of parameter maps. 



Predict fracture paths 



Determine fracture risk 



• Watershed segmentation of selected 



boundaries. 

Search and identify possible fracture paths. 



Evaluate selected parameters along predicted 
fracture paths. 

Generate z-score and t-score for predicted 
fracture paths. 
Estimate fracture risk. 



Figure 18. Process to determine fracture risk using overlapping window processing 
and fracture paths prediction. 
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Figure 22 
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Figure 24 
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